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ABSTRACT 


This  thesis  explores  the  application  of  the  Weibull 
distribution  to  the  Cost  Effectiveness  Analysis 
Spreadsheet  Model  (CEAMOD)  for  the  Aircraft  Engine 
Component  Improvement  Program  (CIP) .  The  current  model 
assumes  the  exponential  distribution  for  all  unscheduled 
failures .  This  thesis  explores  whether  the  Weibull 
distribution  can  be  used  in  the  CEAMOD.  This  process  is 
pg2rformed  through  the  use  of  a  simulation  program  and 
spreadsheet  analysis .  The  thesis  first  examines  past 
concerns  that  the  use  of  the  Weibull  will  create  cyclic 
annual  numbers  of  failures  of  a  component  over  an  engine's 
life.  It  then  examines  the  empirical  distributions  of  the 
number  of  failures  per  year  generated  from  simulations  for 
a  range  of  Weibull  parameters  and  compares  these  with  the 
Poisson  distribution  resulting  from  the  current  CEAMOD 
assumption  of  exponential  times  between  failures . 
Finally,  it  considers  how  life  cycle  costs  considered  by 
the  CEAMOD  will  change  if  the  Weibull  distribution  is 
assumed.  The  conclusion  is  that  the  nature  of  the 
distribution  of  the  annual  number  of  failures  should  be 
studied  further  before  the  Weibull  is  incorporated  into 


the  CEAMOD. 
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I .  INTRODUCTION 


A.  BACKGROUND 

This  thesis  analyzes  the  input  and  methodology  of  the 
statistical  distributions  applied  to  the  General  Electric 
Aircraft  Engines  Cost  Effectiveness  Analysis  Spreadsheet  Model 
(GE  CEAMOD) . 

The  GE  CEAMOD  is  a  spreadsheet  program,  originally 
developed  by  Pratt  and  Whitney,  and  later  refined  by  General 
Electric  for  use  in  analyzing  the  cost  effectiveness  of 
aircraft  engine  change  proposals  [Ref.  1],  The  current 
release.  Version  2.1,  is  written  for  use  with  the  Microsoft 
EXCEL  4.0  program.  Both  the  Navy  and  the  Air  Force  use  the 
CEAMOD  program  in  conjunction  with  their  prime  contractors  to 
determine  whether  an  Engineering  Change  Proposal,  as  submitted 
by  the  manufacturer,  is  cost  effective  over  the  proposed  life 
of  the  engine/component . 

The  CEAMOD  program  is  used  in  support  of  the  aircraft 
engine  Component  Improvement  Program  (CIP) .  The  CIP  has  the 
following  three  objectives  [Ref.  1] : 

1.  To  maintain  an  engine  design  which  allows  the  maximum 
aircraft  availability  at  the  lowest  cost  to  the 
government . 

2.  To  correct,  as  rapidly  as  possible,  any  design 
inadequacy  which  adversely  affects  safety  of  flight. 

3 .  To  correct  any  design  inadequacy  which  causes 
unsatisfactory  engine  operation  or  adversely  affects 
maintainability  and  logistic  support  in  service. 
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The  CEAMOD  program  supports  the  CIP  by  allowing  engineers 
and  program  managers  to  analyze  the  expected  life  cycle  costs 
differential  between  a  current  component  and  a  proposed  change 
in  that  component.  Because  the  life  of  a  component  is 
typically  many  years,  all  future  expected  costs  or  savings  are 
expressed  in  net  present  value  terms .  Since  many  of  the 
factors  included  in  the  analysis  are  estimates,  the  final 
program  output  is  subject  to  interpretation.  In  addition, 
there  are  other  "softer"  factors  that  are  considered  as  well. 
These  include:  safety,  increased  mission  effectiveness,  and 
decreased  overall  downtime  for  maintenance  and  repairs.  The 
CEAMOD  does  not  analyze  these  factors,  but  none-the-less,  they 
are  critical  to  the  overall  decision  of  whether  to  incorporate 
the  modification. 

One  factor  that  has  been  under  discussion  with  the 
program  developers  and  users  is  the  assumption  in  CEAMOD  of  a 
constant  component  failure  rate .  By  assuming  that  the  failure 
rate  of  a  component  is  constant,  the  time  between  failures  can 
be  described  by  the  exponential  probability  distribution. 
There  is  a  belief  that  assuming  this  distribution  ensures  the 
CEAMOD  results,  which  show  fairly  constant  numbers  of  failures 
per  year,  are  representative  of  reality.  However,  the 
designers  and  users  of  CEAMOD  have  proposed  that  the  Weibull 
distribution  may  reflect  the  reality  of  the  times  between 
component  failures  more  accurately.  The  concern  is  that  due 
to  the  Weibull 's  ability  to  reflect  changing  failure  rates. 
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the  CEAMOD  program  may  become  too  complex  [Ref.  9].  In 
addition,  the  question  remaining  to  be  investigated  is:  "can 
the  Weibull  be  successfully  introduced  into  the  CEAMOD  program 
without  creating  a  cyclic  nature  to  the  number  of  component 
failures  over  time?" 

B .  OBJECTIVE 

The  objective  of  this  thesis  is  to  analyze  one  of  the 
specific  assumptions  of  the  CEAMOD  program.  This  assumption 
has  to  do  with  the  current  and  proposed  unscheduled  failure 
rates  for  the  component  in  question.  The  current  program 
assumes  that  the  time  between  failures  of  a  component  is 
distributed  according  to  the  exponential  probability 
distribution.  If  the  time  between  failures  behaves  according 
to  the  Weibull  distribution,  how  are  the  number  of  failures 
per  year  affected?  There  are  two  relevant  issues  with  regards 
to  the  objective,  they  are: 

1.  Can  the  number  of  failures  per  year  be  described  by 
a  simple  probability  distribution? 

2.  Are  there  any  significant  differences  in  the  overall 
life-cycle  costs  for  a  component  when  the 
exponential  distribution  is  replaced  by  the  Weibull? 


C.  SCOPE 

The  scope  of  this  thesis  is  limited  to  an  analysis  of  a 
single  component's  change  and  how  the  use  of  the  Weibull 
distribution  instead  of  the  exponential  affects  the  expected 
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life  cycle  costs.  This  thesis  will  only  focus  on  the  expected 
costs  of  unscheduled  maintenance  (i.e.  unexpected  failures) 
and  not  on  the  costs  associated  with  planned  maintenance.  The 
trial  data  used  was  provided  by  General  Electric  and  was  the 
example  previously  analyzed  in  Rau's  thesis  [Ref.  4] . 


D .  METHODOLOGY 

Due  to  the  nature  and  limitations  of  the  CEAMOD  program, 
a  simple  substitution  of  the  number  of  failures  generated 
assuming  the  Weibull  distribution  was  not  possible  [Ref.  9] . 
The  exponential  distribution  provides  for  a  constant  failure 
rate.  Thus,  if  a  component  reaches  a  certain  point  in  its  life 
where  its  failure  rate  is  constant,  it  has  no  greater  or 
lesser  chance  of  failing  in  the  next  instant  of  time  for  any 
time  during  its  life.  The  number  of  failures  per  year  then  has 
the  Poisson  probability  distribution  with  the  constant  failure 
rate  as  its  mean  number  per  year.  The  Weibull  distribution 
allows  for  a  changing  failure  rate.  The  number  of  failures 
per  year  is  no  longer  Poisson  distributed.  Thus,  before  the 
Weibull  can  be  incorporated  into  the  CEAMOD,  the  nature  of  the 
distribution  of  the  number  of  failures  per  year  must  be 
understood. 

To  study  the  distribution  of  the  number  of  failures  per 
year,  assuming  Weibull  times  between  failures,  a  simulation 
program  was  written  in  the  SIMAN  simulation  language.  A 
version  of  the  CEAMOD  was  also  developed  to  accept  the  SIMAN 
data  so  that  a  financial  analysis  of  the  effects  from  the 
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Weibull  distribution  could  be  conducted. 


E.  ORGANIZATION  OF  THE  THESIS 

Chapter  II  provides  a  background  of  the  current  CEAMOD 
progratn  and  analyzes  its  methodology  and  assumptions .  Chapter 
III  is  an  in  depth  explanation  of  the  proposed  changes  and 
their  expected  effects  on  the  program.  Chapter  IV  describes 
the  development  of  the  SIMAN  program  and  the  financial 
spreadsheet.  Chapter  V  is  an  analysis  of  the  model  output. 
Chapter  VI  contains  a  summary  of  the  research  efforts  and  its 
conclusions . 
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II.  BACKGROUND  OF  THE  CEA  MODEL 


A.  DEVELOPMENT  OF  THE  CURRENT  MODEL 

The  GE  CEAMOD  program  is  an  EXCEL-based  spreadsheet 
program  currently  used  by  the  Navy  and  Air  Force  to  estimate 
the  future  costs  of  aircraft  engine  component  modifications. 
This  program  was  originally  developed  as  a  mainframe  program 
by  Pratt  and  Whitney  at  the  request  of  the  Air  Force.  General 
Electric  later  converted  it  to  a  PC  based  Lotus  1-2-3 
spreadsheet  and  then  to  its  current  configuration  as  an  EXCEL 
spreadsheet  [Ref.  l] . 

The  Navy  and  Air  Force  have  adopted  the  current  model  as 
the  primary  analysis  tool  in  determining  whether  a 
modification  to  an  engine  or  group  of  engines  is  cost 
effective  from  a  life-cycle  point  of  view.  As  noted  in  the 
Reeves  thesis  [Ref.  1],  the  model  does  have  its  imperfections 
and  "illogical"  steps,  but  it  is  currently  the  most  effective 
model  to  use  in  determining  long  term  cost  effectiveness.  The 
CEAMOD  Users  Group  meets  twice  a  year  to  discuss  and  approve 
improvements  to  the  model.  Version  2.2  was  approved  at  the 
meeting  held  on  29  November  1994. 

B.  ASSUMPTIONS  OF  THE  CEAMOD  PROGRAM 

There  are  a  number  of  assumptions  that  have  been 
incorporated  into  the  current  CEAMOD  program.  The  primary 
assumption  is  that  it  will  incorporate  only  one  hardware 
change  at  a  time.  The  significance  of  this  assumption  is  that 
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a  proposed  component  modification  is  considered  to  be  totally 
independent  of  any  other  component.  Thus,  one  modification 
has  no  bearing  on  any  other.  Although  this  may  be  true  in 
Q0]frain  cases,  it  is  not  necessarily  true  in  all  cases.  An 
aircraft  engine  is  a  complex  assembly  and  the  reliability  of 
one  component  can  and  does  have  an  impact  on  other  components . 
This  assumption  was  made  primarily  due  to  cost  considerations 
in  the  development  of  the  CEAMOD  program. 

The  costs  usfed  in  the  model  to  calculate  the  Net  Present 
value  (NPV)  are  primarily  maintenance  and  operational  costs. 
Also  included  in  the  model  are  the  "kit"  costs  and  labor 
required  to  make  the  change.  When  a  component  modification  is 
proposed  and  the  data  is  analyzed  using  the  CEAMOD,  the 
expected  life- cycle  costs  associated  with  the  current 
configuration  are  compared  to  those  of  the  proposed  changes. 
Figure  2.1  is  a  chart  outlining  which  costs  are  included  and 
which  are  not . 

Based  upon  these  assumptions,  the  model  is  designed  to 
ONLY  include  those  logistics  costs  directly  related  to 
maintenance,  operations  and  kit  installation.  All  of  the 
costs  of  developing  the  modification  and  program  management 
are  disregarded.  For  the  most  part  it  is  a  logical  conclusion 
to  not  include  the  sunk  costs  of  R&D  and  D&T,  but  if  there  is 
a  significant  life  cycle  cost  differential  between  the  current 
configuration  and  the  proposed  configuration  for  engineering 
data  and  program  management  these  factors  should  be  included. 
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COST  CATEGORY 


COSTS  INCLUDED  IN  MODEL 


Component  Modification  Cost 


Post-Modification  Material 
Costs 


Cost  of  materials  for 
maintenance  AFTER  modification 


Post-Modification  Labor  Costs  Cost  of  labor  for  maintenance 

AFTER  modification 


Cost  of  re-writing  and 
implementing  new  documentation 


Cost  of  modifying  or  creating 
new  support  equipment 


Cost  of  purchasing  new 
inventory 


Technical  Documentation 
Modification  costs 


Post-Modification  Support 
Equipment  Costs 


Post-Modification  Inventory 
Costs 


Post-Modification  Fuel  Costs  Cost  of  change  in  fuel  usage 


Cost  savings  due  to  prevention 
of  catastrophic  aircraft  loss 


Post-Modification  Aircraft 
Loss  Costs 


COSTS  NOT  INCLUDED  IN  MODEL 


Research  and  Development  Costs  Cost  to  develop  modification 


Funds  spent  to  ensure  that  the 


Design  and  Test  Costs 


Engineering  Data  Costs 


Program  Management  Costs 


component  performs  as  designed 


Cost  to  develop  eng  stds  for 
new  components 


Administrative  costs  related 


to  the  program  management 

Table  2.1:  CEAMOD  cost  Indus ion/exclus ion  table 
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C.  LIMITATIONS  OF  THE  MODEL 

Although  the  CEAMOD  program  is  quite  complex,  there  are 
some  inherent  limitations.  One  critical  limitation  is  the 
spreadsheet ' s  ability  to  incorporate  only  one  component  change 
at  a  time.  This  limitation  is  primarily  due  to  the  program's 
ability  to  accept  only  one  failure  rate.  This  may  be 
overcome,  at  least  partially,  by  combining  several  component 
changes  into  one  "complex"  change.  This  process  of  combining 
failure  rates  is  unrelated  to  the  use  of  the  program,  and  must 
be  calculated  offline  [Ref.  1]  .  Another  possible  way  to 
address  the  problem  is  to  incorporate  one  component  at  a  time 
in  successive  runs  of  the  model.  The  iterative  process  might 
be  to  shift  the  modified  configuration  data  to  the  "current" 
column  and  the  next  modification  to  the  "proposed"  column.  By 
doing  this  for  all  components,  the  overall  cost  of  all 
modifications  can  be  estimated.  This  process  of  incorporation 
must  be  handled  with  care,  because  the  order  of  incorporation 
may  affect  the  outcome.  Some  orderings  may  increase  overall 
costs  because  of  interactions  with  other  parts  of  the  engine. 

Another  limitation  is  related  to  the  input  data  developed 
by  the  contract  or /program  manager.  The  data  that  is  input 
into  the  program,  both  for  costs  and  component  reliability 
must  be  developed  offline.  Much  of  the  data  is  estimated  in 
the  early  stages  of  the  program  rather  than  being  known 
precisely.  Experience  has  shown  that  reality  is  different  than 
the  estimate.  Thus,  this  program  should  not  be  used  only  at 
the  inception  of  the  program,  but  should  be  run  at  any  time 
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there  are  significant  changes  to  the  relevant  data.  If  the 
actual  failure  data  of  new  components  in  service  does  not  meet 
estimates,  there  may  be  a  need  to  re-evaluate  the  modification 
policies.  As  a  program  develops,  estimates  of  reliability 
and  costs  should  be  continually  updated  and  entered  into  the 
program.  Even  small  changes  to  reliability,  cost  or  usage 
estimates  may  have  a  large  effect  on  the  total  life  cycle 
costs.  These  changes  may  lead  the  program  manager  to  search 
for  more  reliable  or  cost-effective  components. 

A  significant  limitation  of  this  program  is  the 
assumption  of  the  exponential  probability  distribution  for 
times  between  failures.  It  will  be  the  focus  of  this  thesis. 
The  assumption  of  the  exponential  distribution  results  in 
expected  failure  rates  being  constant  [Ref.  4].  The 
consequences  of  assuming  the  exponential  distribution  for 
times  between  failure  is  that  the  number  of  failures  over  a 
specified  period  of  time  is  Poisson  distributed  with  a  mean 
which  is  the  product  of  the  failure  rate  and  the  specified 


length  of  time. 

The  input  menu  of  the  CEAMOD  program  requires  a  specific 
number  to  reflect  the  expected  failure  rate  for  both  the 
current  configuration  and  the  proposed  change.  The  rate  is 
expressed  as  the  expected  number  of  failures  per  1000  flight 
hours.  The  definition  of  "failure"  is  intentionally  kept 
vague,  and  must  be  determined  by  the  contractor /program 
manager.  In  the  operating  environment  a  failure  may  be  defined 
as  the  time  at  which  the  component  will  no  longer  function,  or 
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may  be  defined  as  the  point  at  which  it  has  failed  to  meet 
certain  operating  parameters . 

The  assumption  of  a  constant  failure  rate  has  been  a 
topic  of  discussion  for  several  years  by  the  CEAMOD  Users 
Group.  The  reason  is  that  the  reality  of  component  failures 
is  that  they  tend  to  decrease  as  the  "burn-in"  period  passes, 
stabilize  during  the  majority  of  the  expected  life,  and 
increase  in  the  later  "wear-out"  period.  This  phenomenon  can 
be  illustrated  by  the  well-known  "bathtub  curve"  [Ref.  10]  . 
Figure  2.2  is  an  example  of  a  failure  rate  bathtub  curve.  In 
those  periods  of  decreasing  and  increasing  failure  rates,  the 
time  between  failures  are  best  estimated  by  Weibull 
distribution  [Ref.  5].  How  to  incorporate  the  Weibull 
assumption  into  the  CEA  model  has  been  a  question  which  the 
users '  group  has  puzzled  over  for  several  years .  The  purpose 
of  this  thesis  is  to  attempt  to  determine  how  to  do  that 
incorporation.  The  specifics  of  the  Weibull  distribution 
will  be  discussed  in  Chapter  III. 
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RVE 


Figure  2.2;  A  failure  rate  bathtub  curve 
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III.  PROPOSED  CHANGES 


A.  THE  WEIBULL  DISTRIBUTION 

The  Weibull  Distribution  is  named  after  a  Swedish 
Physicist  named  Waloddi  Weibull,  who  in  1939,  developed  it  to 
represent  the  breaking  strength  of  materials  [Ref.  5].  Since 
then  it  has  become  widely  used  in  many  areas  of  analysis, 
particularly  those  related  to  the  modeling  of  reliability  of 
electrical  and  mechanical  components  [Ref.  5].  The  Weibull 
distribution  most  frequently  provides  the  best  fit  for  the 
type  of  failure  data  experienced  in  the  gas  turbine  industry 
[Ref.  5]  .  Thus,  the  consideration  of  using  the  Weibull 
distribution  for  times  between  failures  in  the  CEA  model  has 
been  an  issue  for  the  last  several  years . 

The  probability  density  function  (PDF),  f(t, -11,13),  of  the 
Weibull  distribution  is  represented  by  the  following  equation 
for  values  of  t^O,  ri>0,  and  |3>0  [Ref.  8]  : 


t  p 

r(t;n,P)=4t^'"e 

n 


(3.11 


The  specifics  of  the  Weibull  distribution  are  explained 
below,  preceded  by  definitions  of  its  parameters: 

1.  BetaO)  :  Beta  is  commonly  referred  to  as  the  "slope" 
parameter.  The  value  of  Beta  O>0)  determines  which  member  of 
the  family  of  Weibull  curves  that  best  fits  the  data  being 
analyzed.  The  type  of  failure  may  be  any  of  those  represented 
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by  the  bathtub  curve.  "Slopes"  or  (3  values  less  than  one 
correspond  to  decreasing  failure  rates,  typical  of  those  found 
in  the  infant  mortality  segment  of  the  bathtub  curve, 
"slopes"  or  (3  values  equal  to  one  represent  "random"  failure 
where  the  failure  rate  is  constant,  and  "slopes"  or  p  values 
greater  than  one  correspond  to  increasing  failure  rates, 
typical  of  those  found  in  the  wear-out  segment  of  the  bathtub 
curve.  A  Weibull  distribution  with  a  beta  of  one  is  the  same 
as  the  exponential  distribution. 

2.  Eta(ri):  Eta  defines  the  "characteristic  life"  of  the 
Weibull  distribution.  The  characteristic  life  of  a  component, 
as  represented  by  eta,  represents  the  value  at  which  63.2%  of 
all  failures  will  occur  regardless  of  the  value  of  beta.  The 
value  of  l/h  could  be  defined  as  the  failure  rate  at  the 
characteristic  life  in  the  Weibull  distribution.  The  value  of 
eta  is  the  mean  for  the  Weibull  distribution  (and,  of  course 
the  exponential)  when  the  value  of  beta  is  one,  but  is  not  the 
tnean  when  the  value  of  beta  is  less  than  or  greater  than  one. 
Figure  3.2  shows  graphically  the  shape  of  the  PDF  for  several 
values  of  beta. 

The  hazard  function  describes  the  instantaneous  failure 
rate .  The  hazard  function  for  the  Weibull  distribution  is 
represented  by  the  following  equation  for  value  of  t^O,  h^O 
and  p^O  [Ref .  8]  : 


]3(t;ri,P)  =  {^)^pt^''. 


(3.2) 
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Figure  3.2:  Weibull  PDF  curves 


When  (3  =  1  we  have  the  exponential  function,  the  hazard 

function  becomes: 


Thus,  because  the  failure  rate  equals  l/n,  and  n  equals  the 
mean  for  an  exponential  distribution,  the  hazard  rate  is 
simply  defined  as  one  over  the  mean. 
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The  mean  of  the  Weibull  distribution  for  values  of  beta 
less  than  and  greater  than  one  is  given  by  [Ref .7] : 


P  = 


ri  1 


(3.4) 


where  r(.)  is  the  Gamma  Function.  The  variance  of  the  Weibull 
distribution  is  given  by  [Ref.  7]: 


2 


(3.5) 


A  comparative  example  showing  the  mean  and  variance  for 
different  values  of  beta  is  provided  in  table  3.1: 


p 

h 

t 

h(t) 

MTBF 

2 

.5 

50 

50 

.01 

100 

50000 

1 

50 

50 

.02 

50 

2500 

2 

50 

50 

.04 

44.3 

536.5 

Table  3.1:  Comparison  of  means  and  variances 


Table  3 . 1  demonstrates  that  even  though  each  example  has 
the  same  t\  value  (characteristic  life)  the  actual  mean  values 
are  very  different.  With  a  (3  =  1,  or  an  exponential 
distribution,  the  mean  is  equal  to  the  h  value.  But,  when  the 
beta  value  is  less  than  one  the  mean  is  much  higher  and  when 
the  beta  value  is  greater  than  one,  the  mean  is  less  than  the 
p  value.  The  variance  also  decreases  as  the  value  of  beta 
increases . 
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B.  DETERMINATION  OF  WEIBULL  PARAMETERS 

In  order  to  properly  determine  the  Weibull  distribution 
which  fits  a  given  set  of  failure  data,  the  values  of  beta  and 
eta  need  to  be  estimated.  This  process  is  called  Weibull 
analysis  [Ref.  5]  . 

Weibull  analysis  is  best  performed  using  Weibull  graph 
paper,  an  example  of  which  is  in  Appendix  A.  In  performing  a 
Weibull  analysis,  there  are  several  questions  that  need  to  be 
answered.  The  primary  one  being,  "Can  the  data  be  described 
by  the  Weibull  distribution?"  A  proper  Weibull  plot  is 
represented  by  a  straight  line,  thus  if  the  data  does  create 
a  straight  or  approximately  straight  line  (where  the  majority 
of  the  points  fall  in  a  straight  line,  but  there  may  be  a  few 
outliers)  ,  then  the  data  can  be  represented  by  the  Weibull 
distribution.  Appendix  A  also  contains  a  typical  plot  of 
Weibull  data  on  Weibull  graph  paper. 

Weibull  graph  paper  has  unique  scales  for  the  X  and  Y 
axes,  which  graphically  will  create  a  straight  line  if  the 
data  is  Weibull.  The  values  on  the  X  axis  represent  the 
measure  of  life  for  the  data,  and  are  used  as  the  primary 
estimator  for  the  value  of  eta.  The  Y  axis  represents  the 
expected  cumulative  percentage  of  components  that  have  failed. 
By  measuring  from  the  63.2%  line  on  the  Y  axis,  the  value  of 
beta  can  be  determined. 

Once  time  to  failure  data  has  been  collected,  it  must  be 
ordered  from  first  failure  to  last  failure.  Once  in  the 
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propsr  order,  the  Median  rank  inust  be  determined.  Appendix  B 
provides  an  example  of  a  table  for  the  determination  of  the 
median  rank  based  upon  the  rank  order  and  the  sample  size. 
Appendix  B  depicts  the  following  (Sample  size  =  3) : 

Time  to  Failure  Median  Rank 

1  =  100  hrs  20.6 

2  =  200  hrs  50.0 

3  =  300  hrs  79. 3 

These  values  are  plotted  on  the  Weibull  paper,  with  the 
hours  on  the  X  axis  and  the  median  rank  on  the  Y  axis.  Next, 
a  straight  line  is  fitted  to  the  data.  Then  eta  is  determined 
by  drawing  a  vertical  line  down  to  the  X  axis  from  the  point 
where  the  data  line  crosses  the  Y=63.2%  horizontal  line.  This 
point  on  the  X  axis  is  the  value  of  the  characteristic  life. 
In  the  case  of  the  chart  in  Appendix  A,  that  value  is  roughly 
240  hours. 

Beta  is  determined  by  moving  down  the  data  line  to  a 
point  which  is  horizontally  one  inch  away  from  the  vertical 
line  corresponding  to  the  characteristic  life  and  then 
measuring  down  from  the  intersection  of  the  63.2%  line  and  the 
data  line  to  the  intersection  of  the  horizontal  and  vertical 
lines.  Slope  is  defined  as  rise  over  run,  thus  the  distance 
in  inches  down  to  the  "one-inch  horizontal  intersection"  is 
the  slope  or  the  Beta  value.  In  Appendix  A  the  distance  is 
1.75  inches  and  therefore  the  value  of  Beta  is  1.75.  Because 
the  value  of  Beta  is  greater  than  one,  this  indicates  that  the 
data  is  from  a  component  which  in  the  wearout  stage  of  the 
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bathtub  curve . 


Figure  3.5  depicts  the  probability  density 
function  of  this  data,  with  a  p  =  1.75  and  an  n  =  240. 


Figure  3.3:  Weibull  PDF  for  Appendix  A  example's 
pareuneters . 

C.  PROPOSED  CHANGES  TO  THE  CEAMOD  PROGRAM 

The  proposal  to  incorporate  the  Weibull  distribution  into 
the  CEAMOD  program  has  been  studied  by  General  Electric  and 
reference  9  states  that  they  feel  that  Weibulls  are  too 
complex  for  this  model.  However,  a  Weibull  distribution  has 
the  capability  of  allowing  a  failure  rate  that  will  change 
with  time.  For  a  given  component  in  a  system,  the  Weibull 
allows  either  an  increasing,  constant,  or  decreasing  failure 
rate.  In  contrast,  the  exponential  distribution  assumes  only 
that  the  failure  rate  is  constant.  Thus,  if  a  component  has 
reached  a  given  number  of  operating  hours,  there  is  no  greater 
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or  lesser  chance  of  failure  after  this  point  in  time  than  at 
any  other  point  in  its  life. 

From  an  analytical  point  of  view  it  is  extremely 
difficult  to  incorporate  Weibull  generated  failures  into  the 
CEAMOD  program.  The  reason  is  that  the  distribution  of 
failures  over  a  given  time  interval  is  difficult  to  derive 
except  in  the  case  of  p  =  1.  For  that  case  the  Poisson 
distribution  can  be  used  since  times  between  failures  are 
exponentially  distributed.  However,  Professor  James  Esary  of 
the  Operations  Research  Department  of  the  Naval  Postgraduate 
School  suggested  that  the  use  of  simulation  results  should 
provide  an  empirical  distribution  which  might  be  useful  in  the 
CEA  model.  Thus,  this  thesis  uses  a  program  to  simulate  what 
might  happen  if  a  component's  times  between  failure  were 
Weibull  distributed. 

This  thesis  will  analyze  simulations  of  a  component's 
failures  over  a  period  of  the  expected  life  of  an  engine.  To 
provide  a  basis  for  comparison,  the  general  operating 
parameters  for  this  component  are  assumed  to  be  those  provided 
by  the  General  Electric  example  (a  copy  of  the  CEAMOD  output 
for  that  example  is  presented  in  Appendix  C) .  One  intent  of 
comparing  the  behavior  of  component  failures  represented  by 
the  exponential  distribution  to  those  represented  by  the 
Weibull  distribution,  is  to  determine  whether  there  is  a 
significant  cost  difference,  in  the  long  term,  based  upon  the 
net  present  value  of  expenditures  and  savings . 
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IV.  SIMULATION  AND  SPREADSHEET  MODEL  DEVELOPMENT 


A.  MODEL  CONCEPT 

The  concept  for  the  model  is  to  simulate  failures  of  a 
component  over  its  expected  life.  The  model  incorporates 
failures  that  occur  either  in  an  exponential  manner  or  in  a 
Weibull  manner  over  a  period  of  ten  years.  As  noted  in  the 
previous  chapter  by  determining  the  time  to  failure  and  the 
annual  number  of  failures,  it  should  be  possible  to  determine 
whether  the  use  of  the  Weibull  distribution  is  feasible. 


B.  SIMAN  SIMULATION  PROGRAM 

In  order  to  simulate  this  process,  a  simulation  program, 
written  in  the  SIMAN  language  was  developed.  Appendix  D 
contains  the  program  used  to  create  the  failure  data. 

The  program  begins  by  introducing  ten  engines  into  the 
system  at  time  zero.  It  is  assumed  that  all  ten  engines  are 
in  a  "ready  for  flight"  condition.  The  component  of  interest 
fails  based  upon  the  assumed  distribution.  The  engine  then 
enters  into  a  repair  facility,  is  repaired,  and  the  engine  is 
then  again  ready  for  service.  That  engine  fails  again  at  a 
future  time  based  upon  the  assumed  distribution,  and  the 
process  is  repeated. 

In  order  to  run  this  program,  several  assumptions  were 
made.  These  are  shown  in  Table  4.1.  The  basic  time 
parameters  were  taken  from  the  data  provided  by  the  General 
Electric  example,  and  others  were  assumed  by  the  author. 
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PARAMETER 

VALUE 

■MfUif 

ANNUAL  FLIGHT  HRS 

240 

GE  Example 

MEAN  TIMES  TO  FAIL 

20,50,100,500 

GE  Example 

(HOURS) 

Author 

REP  FAC  DELAY  TIME 

0 

Author 

(HOURS) 

NBR  ENGINES /YEAR 

10 

Author 

NBR  OF  REP  FAC 

1 

Author 

LIFE  OF  ENG  (YRS) 

10 

Author 

Table  4.1:  SIMAN  parameters  for  the  failure  event  simulation. 

The  General  Electric  example  uses  an  average  annual 
fight  hours  value  of  240  hours  per  aircraft  per  year.  For  the 
component  being  analyzed  in  the  example,  the  assumption  is 
that  the  current  component  has  a  MTBF  of  50  hours  and  the 
proposed  component  has  a  MTBF  of  500  hours.  In  addition  to 
these  values,  MTBF  values  for  the  exponential  distribution  of 
20  and  100  were  chosen  as  well.  The  reason  for  these 
additional  values  was  to  simulate  potential  failures  at  a 
higher  failure  rate  (20)  ,  and  those  at  a  rate  (100)  between 
the  two  values  of  the  example  problem.  The  benefit  of  this  is 
to  provide  additional  information  for  failure  rate 
characteristics  not  necessarily  in  line  with  the  original 
assumptions.  The  eta  value  (characteristic  life)  for  the 
Weibull  distribution  was  assigned  using  the  same  values  as  for 
the  mean  of  the  exponential  distribution.  The  beta  factors 
were  all  greater  than  one  corresponding  to  a  component  having 
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an  increasing  failure  rate.  The  beta  values  ranged  from  1.5 
to  5  (see  Table  4.2) . 

The  repair  facility  delay  time  of  zero  hours  was  a 
necessary  assumption  to  ensure  the  SIMAN  program  was 
comparable  with  the  CEAMOD  in  its  assumptions  of  repair  time. 
The  number  of  engines  introduced  per  year,  ten,  and  the  engine 
life  value,  ten  years,  were  chosen  to  provide  sufficient  data 
points  to  complete  a  credible  analysis .  The  number  of  repair 
facilities,  one,  has  no  impact  on  this  simulation,  but  if  a 
repair  delay  time  was  introduced  it  could  impact  the  number  of 
failures  per  year  by  being  a  potential  bottleneck. 

Once  the  parameter  selection  was  made  and  the  program 
written,  the  probability  distribution  of  times  between 
failures  was  assigned.  As  discussed  above  the  author  chose 
the  following  exponential  and  Weibull  distributions  to  test  in 
the  program: 


DISTRIBUTION 

EXPO(MTBF),  WEIBULL(ri,  3) 

EXPONENTIAL 

20,50,100,500  (hours) 

WEIBULL 

(20,1.5)  (20,2)  (20,3)  (20,5) 

WEIBULL 

(50,1.5)  (50,2)  (50,3)  (50,5) 

WEIBULL 

(100,1.5) 

(100,2) 

(100,3) 

(100,5) 

WEIBULL 

(500,1.5) 

(500,2) 

(500,3) 

(500,5) 

Table  4.2:  Distributions  and  their  parameters  used  for  failure 
times  in  the  simulation  (all  times  are  expressed  in  hours) . 
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As  discussed  in  Chapter  III,  the  actual  means  and 
variances  for  the  Weibull  are  different  than  those  for  the 
exponential  for  the  same  eta  value .  Table  5 . 3  reflects  the 
theoretical  means  and  variances  in  hours  for  the  assumed 
distributions  presented  in  table  4.2.  These  values  for  the 
means  and  variances  were  computed  using  the  equations  (3.4) 
and  (3.5)  from  Chapter  III.  The  theoretical  means,  for  any 
value  of  eta,  tend  to  decrease  as  beta  increases  until  it 
reaches  approximately  2.2.  At  this  point  the  mean  begins  to 
increase  slowly,  asymptotically  approaching  the  mean  of  the 
exponential.  The  variances,  for  all  values  of  eta,  decrease 
as  beta  increases.  Theoretically,  for  a  beta  value  of 
infinity,  the  mean  for  the  Weibull  equals  the  mean  for  the 
exponential  and  the  variance  is  zero. 

The  SIMAN  program  begins  by  creating  ten  "entities" 
(engines)  and  assigns  consecutive  engine  numbers,  "ENGNUMs", 
to  each,  beginning  at  one.  The  time  of  arrival  into  the 
program  is  assigned  as  the  current  time  in  the  simulation. 
With  the  ten  engines,  all  operationally  capable,  the  program 
"delays"  each  engine  by  a  time  equivalent  generated  by  the 
assumed  distribution.  This  delay  simulates  the  operational 
phase  of  each  engine.  At  the  end  of  the  delay  the  engine 
component  fails. 

Once  a  failure  has  occurred,  the  engine  enters  into  a 
repair  queue  to  await  repair.  Once  in  repair  the  author 
assumes  that  the  time  in  the  repair  facility  is  zero. 
Although  this  is  not  a  realistic  assumption  for  real  world 
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operations,  it  was  necessary  to  approximate  the  assumptions  of 
the  CEAMOD  program.  The  CEAMOD  assumes  that  each  engine  fails 
according  to  the  exponential  distribution,  and  does  not  assume 
that  the  repair  channel  may  be  a  delaying  factor  in  the  return 
to  operational  capability.  Thus  the  author  has  chosen  a  delay 
of  zero  for  the  repair  channel.  This  SIMAN  program  can,  of 
course,  be  easily  modified  to  reflect  the  reality  of  repair 
channel  delays . 

After  the  repaired  engine  exits  the  repair  channel,  the 
time  to  failure  is  tallied  for  each  engine.  The  engine  is 
then  returned  to  the  point  where  the  next  delay  time  is 
generated.  It  is  then  delayed  again  by  that  amount  and  the 
process  repeats. 

Based  upon  the  assumptions  stated  in  Table  4.1,  the 
program  runs  for  a  total  of  2400  hours  for  each  engine;  that 
is,  240  hours  per  year  for  ten  years.  Because  there  is  no 
"delay"  for  each  engine  in  the  repair  channel,  the  "system"  is 
in  a  steady  state  with  no  need  for  an  initial  warm-up  period. 
If  a  delay  time  for  repair  were  introduced,  the  SIMAN  program 
would  need  to  include  some  time  for  a  warm-up  in  order  to 
achieve  steady  state. 

The  data  that  was  written  to  the  tally  files  reflects  the 
point  in  time  when  each  failure  occurred.  An  example  output 
from  the  program  is  included  in  Appendix  E.  This  output 
provides  the  program  user  with  several  important  pieces  of 
data.  For  each  engine  it  lists  the  average  time  between 
failures  (MTBF) ,  the  variation  of  the  simulated  times,  the 
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Distribution 


MTBF  (Hours) 


Variance  (Hours) 


Expo (20) 


20.0 


400.0 


Weib (20,1.5! 


18 . 1 


150.3 


Weib (20 ,2) 


17.7 


85.8 


Weib (20, 3; 


17 . 9 


42.1 


Weib (20 , 5) 


18.4 


17.7 


Expo (50) 


50 . 0 


2500 . 0 


Weib(50, 1.5; 


Weib (50,2! 


45 . 1 


44.3 


939.2 


536.5 


Weib (50, 3) 


44.6 


263.3 


Weib  (50,5) 


45.9 


110.6 


Expo (100] 


Weib (100,1.5) 


100.0 


90.3 


10000.0 


3756.9 


Weib (100,2! 


88 . 6 


2146 . 0 


Weib (100, 3) 


89 . 3 


1053.3 


WeibdOO,  5) 


Expo (500) 


91.4 


500.0 


442 . 3 


250000 . 0 


Weib(500,1.5) 


451.4 


93922 . 6 


Weib (500,2) 


443.1 


53650.5 


Weib(500, 3) 


Weib(500,5) 


446.5 


459.1 


26333.2 


11057.5 


Table  4.3;  Theoretical  means  and  variances  for  exponential 
and  Weibull  distributions . 
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minimum,  maximum,  and  the  total  numher  of  failures  that  could 
be  expected  over  the  life  of  2400  hours. 

The  program  also  writes  the  actual  failure  times  to  "dat" 
files,  as  specified  in  the  experimental  frame.  These  "dat" 
files  are  important  in  analyzing  the  behavior  of  the  failures, 
and  the  effect  of  the  assumed  statistical  distribution.  An 
example  of  the  output  from  one  file  is  included  in  Appendix  F . 
The  "dat"  files  can  be  converted  to  a  spreadsheet  readable 
format  by  using  the  "export"  feature  of  SIMAN. 

After  importing  the  data  into  a  spreadsheet,  it  can  be 
easily  analyzed  and  graphed  to  reflect  the  potential  number  of 
failures  each  year  over  the  life  of  the  engines.  By  reading 
the  data  from  all  ten  engines  into  the  spreadsheet,  and 
grouping  the  number  of  failures  by  years  (group  together  the 
failures  from  all  engines,  based  upon  the  time  of  failure)  it 
is  easy  to  determine  how  many  engines  are  expected  to  fail  in 
any  given  year .  These  groupings  can  then  be  graphed  to 
further  visualize  the  number  of  failures  each  year  over  the 
life-cycle.  Examples  of  such  graphs  are  included  in  Appendix 
G. 

The  use  of  a  simulation  program  is  one  way  to  generate  a 
distribution  of  failures,  and  not  necessarily  the  easiest  way. 
By  using  a  spreadsheet,  such  as  Excel,  a  distribution  can  be 
created  for  the  Weibull  distribution.  The  benefit  of  this 
method  is  that  it  can  be  easily  accomplished  without  the 
complication  or  expense  of  a  simulation  program.  The  drawback 
is  that  other  parameters,  such  as  a  repair  channel  delay  time 
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cannot  be  introduced. 

The  process  for  creating  the  Weibull  failure  times  can  be 
accomplished  with  the  following  equation:. 

Weibulhx]  . 

The  value  (-In (rand))  is  a  spreadsheet  command  that  takes  the 
negative  of  the  natural  log  of  a  randomly  generated,  uniformly 
distributed  number  between  zero  and  one.  By  copying  this 
equation  down  a  column  in  the  spreadsheet,  an  entire  set  of 
failures  can  be  created.  This  equation  creates  a  time  to 
failure  for  each  cell.  Column  D  of  table  4.4  is  an  example  of 
the  output  from  this  process.  In  that  table  column  E  is  then 
used  to  show  the  cumulative  clock  time.  This  is  necessary  to 
determine  the  cut-off  point  for  failures  over  a  specific 
period  of  time.  Counting  up  the  number  of  failures  from 
Column  D  provides  the  number  of  failures  over  a  given  period 
of  time.  From  Table  4.4  it  is  easy  to  see  that  during  a 
period  of  240  hours,  there  were  seven  failures.  That  is 
determined  by  simply  counting  the  number  of  failures  that 
occurred  prior  to  the  "cut-off"  time  of  240. 

In  order  to  analyze  many  engines  over  many  years,  the 
equation  for  generating  failures  can  be  copied  to  multiple 
columns,  one  for  each  engine .  To  determine  the  distribution 
of  number  of  failures  over  all  engines  per  year,  the  failure 
time  columns  can  be  analyzed  using  the  "data  analysis" 
function  in  Excel,  and  histograms  representing  number  of 
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failures  in  a  year  can  be  plotted.  This  process  is  quite 
easy,  and  by  changing  the  values  of  beta  and  eta,  the 
distributions  for  all  Weibull  parameters  can  be  determined. 


A 

B 

C 

D 

E 

F 

1 

BETA  =  1 

Failure 

Cum  Time 

2 

ETA  =  50 

28.97 

28.97 

3 

19.77 

48.74 

4 

14.33 

63.07 

5 

9.76 

72 . 83 

6 

54.04 

126.87 

7 

32.16 

159.03 

8 

62.51, 

221.54 

9 

58.78 

280.32 

Table  4.4;  Example  of  spreadsheet  generated  failure  times. 


In  order  to  determine  the  potential  effects  of  a  given 
distribution,  relative  to  the  current  CEAMOD  program's 
assumptions,  the  data  must  be  analyzed  from  a  financial 
perspective.  That  is  where  a  modified,  simplified  spreadsheet 
analysis  comes  into  play.  The  next  section  will  address  this 
function. 

C.  EXCEL  FINANCIAL  SPREADSHEET  PROGRAM 

The  financial  spreadsheet  developed  by  the  author 
computes  the  net  present  value  (NPV)  of  the  proposed  changes 
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over  the  life  of  the  program.  Appendix  I  contains  examples 
of  the  spreadsheet  program.  The  following  are  the  specific 
assumptions  and  calculations  included  in  the  program. 

1.  ASSUMPTIONS  AND  CRITICAL  VALXJES  (TOP  LEFT  CORNER  OF 
SPREADSHEET) 

a.  Interest  Rate  (Row  1) 

The  assumed  interest  rate  is  10%  and  is  the  discount 
rate  at  which  the  cost  or  savings  from  the  modification  will 
be  valued  for  all  future  years.  Thus  a  cost  of  $1000.00  next 
year  is  equivalent  to  $909.09  today.  This  is  obtained  by 
dividing  the  cost  of  $1000.00  by  (1  +  interest  rate) .  For  all 
future  years,  (1  +  interest  rate)  is  raised  to  the  power  of 
the  number  of  years.  Thus  the  value  of  $1000.00  in  five  years 
is  equivalent  to  1000/(1.10)®  or  $620.92  today. 

b.  Flight  Hours  per  Year  (Row  2) 

This  value  is  defined  as  the  expected  number  of 
flight  hours  each  engine  will  operate  for  in  each  calendar 
year.  Thus,  for  example,  with  ten  engines  operating,  that 
value  is  multiplied  by  ten. 

c.  Cost  of  Modification  (Row  3) 

This  cost  is  the  expected  cost  of  the  modified 
component  which  is  in  excess  of  the  current  component's 
replacement  cost .  This  value  can  be  positive  or  negative 
depending  upon  the  specific  component.  This  value  does  not 
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includs  r©s6a.3rch  and  devslopitisnt  costs .  It  is  simply  th.6 
(jif f 0]^0];iQ0  p0twssn  ths  cost  of  a  nsw  unit  of  ths  old  componsnt 
and  the  cost  of  a  new  unit  of  the  new  component. 

d.  Labor  Cost  (Row  4) 

The  cost  per  hour  to  have  trained  technicians 
replace  the  failed  component.  This  price  is  an  aggregate  for 

all  technicians  needed  to  perform  the  job.  This  labor  cost 

/ 

may  be  at  either  the  organizational  level,  intermediate  level, 
or  depot  level  or  all  three.  The  assumption  of  $25.00  per 
hour  is  the  organizational  level  labor  rate. 

e.  Current  Failure  Rate  (Row  5) 

This  rate  is  the  expected  number  of  failures  per 
engine  per  flight  hour.  It  is  equivalent  to  one  divided  by 
the  Mean  Time  Between  Failures  in  flight  hours  for  the 
exponential  distribution.  Therefore  if  the  failure  rate  is 
.02  per  flight  hour  then  the  MTBF  is  50  flight  hours.  With 
240  flight  hours  per  year,  the  expected  number  of  failures  is 
(.02  *  240)  or  4.8  failures  per  year  assuming  the  exponential 
distribution  for  times  between  failures .  The  value  displayed 
in  Appendix  I  is  the  expected  number  of  failures  per  flight 
hour,  rounded  to  three  significant  digits. 


f.  Proposed  Failure  Rate  (Row  6) 

The  proposed  failure  rate  is  defined  in  the  same  manner 
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as  the  current  failure  rate.  It  represents  the  failure  rate 
of  the  proposed  design  change.  This  value  is  used  to 
determine  the  total  expected  failures  of  the  proposed 
configuration  per  year,  and  is  rounded  to  three  significant 

digits . 


g.  LaJsor  Hours  (Row  7) 

This  represents  the  number  of  hours  of  labor  which  will 
be  reguired  to  repair  an  engine  once  it  has  failed.  The  labor 
hours  multiplied  by  the  Labor  cost  per  hour  represents  the 
total  cost  (of  labor)  to  repair  each  engine.  For  the  example 
in  Appendix  I,  ten  is  the  assumed  number  of  labor  hours.  This 
value  could  be  incorporated  into  the  SIMAN  model  as  the  repair 
channel  delay  time. 

h.  Current  Failures  per  Year  (Row  8) 

This  value  represents  the  total  number  of  expected 

failures  per  year  per  engine  based  upon  the  current 

configuration.  This  value  is  obtained  from  the  output  of  the 
SIMAN  program  by  taking  the  average  number  of  failures  for  all 
engines  over  the  ten-year  period.  If  there  are  1000  failures 
for  ten  engines  in  ten  years,  that  equates  to  an  average  of 
ten  failures  per  engine  per  year. 

i.  Projected  Failures  per  Year  (Row  9) 

This  value  represents  the  total  number  of  expected 

failures  per  year  per  engine  based  upon  the  proposed 
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QQnf iguifat ion .  This  vslus  is  obtsinsd.  firotn  ths  output  of  tho 
SIMAN  program  and  is  determined  in  the  same  manner  as  the 
current  failures  per  year. 

2 .  Spreadsheet  Calculations 

a.  Year  (Column  1) 

This  is  the  designation  of  the  year  in  the  life-cycle 
analysis:  it  begins  at  year  one  and  continues  through  year 
twenty.  Based  upon  the  assumed  engine  life,  an  engine 
introduced  in  year  one  will  be  phased  out  of  service  ten  years 
later,  in  year  eleven. 

b.  Engines  Introduced  (Column  2) 

This  represents  the  number  of  engines  with  the  improved 
component  introduced  in  a  given  year.  For  sake  of  simplicity 
the  author  has  chosen  to  introduce  ten  engines  in  each  of  the 
first  ten  years  of  the  cycle.  Each  engine  component  is 
assumed  to  be  introduced  at  the  beginning  of  the  year. 

c.  Attrition  (Column  3) 

This  represents  the  number  of  engines  or  engines  with 
specific  components  removed  from  service  in  a  given  year.  The 
author  has  assumed  a  ten-year  engine  life  and  therefore  the 
calculations  show,  for  example,  the  engines  being  introduced 
in  year  one  are  being  removed  at  the  start  of  year  eleven. 


35 


d.  Total  in  Service  (Column  4) 

This  represents  the  total  number  of  engines  operating 
during  a  given  year.  It  is  the  suiri  of  all  previous 
introductions  minus  all  previous  attritions .  All 
introductions  and  attritions  are  assumed  to  occur  at  the  start 
of  the  year. 


e.  Total  Flying  Hours  (Column  5) 

This  represents  the  expected  total  number  of  flying  hours 
performed  in  a  given  year  for  all  operational  aircraft.  This 
number  is  obtained  by  multiplying  the  expected  number  of 
flying  hours  per  aircraft  per  year  by  the  total  number  of 
engines  in  service  during  that  year. 

f.  Current  Expected  N\nnber  of  Failures  (Column  6) 

This  represents  the  total  number  of  failures  that  can  be 
expected  in  a  given  year  for  the  current  configuration.  This 
value  is  obtained  by  multiplying  the  total  flying  hours  value 
by  the  current  failure  rate  (row  8  of  the  assumption  section) . 
It  may  appear  that  the  value  in  column  6  does  not  match  this 
calculation  because  the  value  displayed  for  the  current 
failure  rate  is  only  displayed  to  three  significant  digits. 
The  actual  value  used  is  the  complete  value  of  the  current 
failure  rate  (ten  digits)  without  rounding. 

g.  Proposed  Expected  Number  of  Failures  (Column  7) 

This  value  represents  the  total  number  of  failures  that 
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can  be  expected  in  a  given  year  for  the  proposed 
^on .  This  value  is  obtained  in  the  satue  manner  as 
column  6  except  that  the  "proposed"  value  (row  9  of  the 
assumptions  section)  is  used  vice  the  "current"  value  (row  8 
of  the  assumptions  section) . 

h.  Current  Costs  (Column  8) 

This  value  represents  the  current  year  dollar  value  of 
expected  outlays  for  maintenance  of  the  current  configuration 
due  to  failures.  This  figure  is  obtained  by  multiplying  the 
Current  Expected  Failure  values  (column  6)  by  the  labor  cost 
and  the  labor  hours  value.  In  the  first  example  in  Appendix 
I,  for  year  one,  the  current  cost  equals  (120.7  failures  x  $25 
per  hour  x  10  hours)  or  $30,175.00.  The  sum  at  the  bottom  of 
the  column  represents  the  total  costs  over  the  entire  life 
cycle  of  operating  an  engine  without  a  component  modification. 

i.  Projected  Costs  (Column  9) 

This  value  represents  the  current  year  dollar  value  of 
the  expected  outlays  for  maintenance  due  to  failures  of  the 
proposed  configuration.  It  is  obtained  by  multiplying  the 
number  of  engines  introduced  (column  2)  by  the  cost  of  a 
modification,  and  adding  that  to  the  projected  expected 
failures  (column  7)  multiplied  by  the  product  of  labor  cost 
and  labor  hours.  The  entry  in  the  first  example  in  Appendix 
I  for  year  one  equals 

(10  new  engines  x  $1000  mod  cost)  +  (47.2  expected  failures  x 
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$25.00  labor  cost  x  10  labor  hours)  which  totals  $21,800.00. 
The  sum  at  the  bottom  of  the  column  represents  the  total 
projected  costs  over  the  entire  life  cycle  of  operating  an 
0ngine  with  a  component  modification. 


j .  Projected  Savings  (Column  10) 

This  value  is  the  difference  between  the  Current  Costs 
and  the  Projected  Costs  (column  8  minus  column  9) .  This  value 
may  be  either  positive  or  negative  depending  upon  the  overall 
cost  differential.  The  sum  at  the  bottom  of  this  column 
represents  the  total  savings  over  the  entire  life  cycle  of 
operating  engines  with  the  modification  instead  of  without  the 
modification . 


k.  NPV  of  Projected  Savings  (Column  11) 

This  value  is  the  cumulative  discounted  value  of  the 
projected  savings  in  each  year.  Since  all  costs  are  assumed 
to  occur  at  the  beginning  of  the  year,  the  first  year's  value 
is  not  discounted.  The  present  value  of  second  year's  cost 
represents  the  value  of  the  projected  savings  for  year  two 
divided  by  (1  +  interest  rate)y^^^-\  plus  the  savings  from  the 
previous  year.  In  the  first  example  in  Appendix  I,  the 
savings  for  year  three  equals  (45,125/(1.1)^  +$32,693.00  = 
$69,986.57.  The  value  at  the  bottom  of  the  column  represents 
the  total  net  present  value  of  the  savings  from  making  the 
configuration  change.  This  is  usually  the  key  value  of  the 
CEA  model .  This  value  is  used  by  the  decision  maker  in 
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deciding  whether  it  is  worthwhile  to  make  the  change. 


D .  SUMMARY 

By  using  the  SIMAN  program  developed  by  the  author,  and 
the  EXCEL  spreadsheet  program,  it  is  possible  to  determine 
whether  there  can  be  any  expected  differences  between  the  use 
of  the  Exponential  and  Weibull  distributions. 

The  simulation  of  the  life  cycle  of  engine  failures  is 
critical  to  this  analysis.  By  looking  at  the  simulated 
failures  of  engines  based  upon  an  assumed  distribution,  it  is 
possible  to  project  the  future.  In  an  operational  arena, 
through  the  use  of  Weibull  analysis,  the  engineer  would,  with 
reasonable  certainty,  know  what  the  failure  distribution  would 
be.  By  then  using  this  distribution  in  a  simulation  program, 
he/she  could  see  what  would  be  expected  in  the  long  run. 

The  CEAMOD  program  has  assumed  that  the  distribution  is 
exponential,  even  though  the  Weibull  may  be  a  more  logical 
choice.  Components  of  mechanical  equipment  do  tend  to  wear 
out  as  time  goes  by.  The  Weibull  allows  the  analyst  to 
reasonably  determine  what  is  the  expected  life  remaining  of  a 
component  once  it  has  operated  for  a  given  number  of  hours. 

In  the  next  chapter,  the  output  data  from  the  simulation 
program  and  the  financial  spreadsheet  will  be  analyzed,  and  a 
comparison  will  be  made  considering  the  various  parameter 
changes.  For  example,  by  comparing  the  expected  number  of 
failures  from  the  exponential  and  the  Weibull  distributions, 
over  the  life  cycle  of  the  engines,  it  can  be  determined 
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whether  there  is  any  merit  to  changing  to  a  Weibull -based 
analysis . 
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V.  ANALYSIS  OF  OUTPUT 


A .  INTRODUCTION 

This  chapter  will  focus  on  the  output  provided  by  the 
SIMAN  and  financial  spreadsheet  programs.  The  purpose  of  this 
section  is  to  compare  the  output  data  from  the  simulation 
using  the  CEAMOD  assumption  of  an  exponential  distribution 
versus  a  Weibull  distribution. 

B.  SIMAN  PROGRAM  OUTPUT 

In  order  to  produce  the  necessary  output  data  for  each  of 
the  relevant  distributions  (see  Table  4.2),  it  was  necessary 
to  update  the  program  after  each  run.  The  only  change  which 
needed  to  be  made  was  to  modify  the  "Delay"  factor  to  reflect 
the  distribution  needed  to  produce  the  results. 

The  first  analysis  was  to  look  at  the  simulated 
MTBFs  of  the  components  over  the  expected  life  of  the  engines. 
In  order  to  determine  the  validity  of  the  simulation  model  it 
was  necessary  to  compare  the  SIMAN  results  to  the  theoretical 
means  and  variances.  Table  5.1  provides  that  comparison. 
Based  upon  the  comparison  of  the  means  and  variances  from  a 
simulation  of  100  engine  years  (ten  engines  for  ten  years 
each)  to  the  means  and  variances  of  the  theoretical 
distributions,  it  appears  that  the  simulation  replicated 
theory  very  closely.  The  highest  difference  between  means 
occurred  for  Weib(500,2)  but  the  simulated  mean  was  only  5.7% 
lower  than  the  theoretical  mean.  The  highest  variance 
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E  (20) 


W(20, 1.5) 


W(20, 3) 


W(20, 5) 


E  (50) 


W(50, 1.5) 


W(50, 2) 


W(50, 3) 


W(50, 5) 


E  (100) 


W (100, 1.5) 


W(100, 2) 


W(100, 3) 


W(100, 5) 


E  (500) 


W(500, 2) 
W(500, 3) 
W(500, 5) 


OBS 

MEAN 


19 . 8 


18.0 


17 . 9 


18.4 


47.1 


45.5 


44.3 


44 . 9 


46.1 


108 . 0 


95.3 


92.0 


91.6 


93.3 


517.2 


419 . 1 

447 . 8 

458 . 8 


THEOR 

DIFF 

OBS 

THEOR 

MEAN 

(0-T) 

VAR 

VAR 

20.0 

(  .2) 

355 . 7 

400 

18.1 

.1 

143 . 0 

150.3 

17 . 7 

0 

82.7 

85.8 

17 . 9 

0 

40.4 

42.1 

18.4 

0 

16.6 

17 . 7 

50 . 0 

(2.9) 

1964.4 

2500 

45 . 1 

.4 

890.9 

939.2 

44.3 

0 

521.6 

536.5 

44 . 6 

.3 

259 . 8 

263.3 

45 . 9 

.2 

108.3 

110.6 

100.0 

8 

9967.3 

10000 

90.3 

5 

3658.3 

3756 . 9 

88.6 

3.4 

2056.3 

2146 . 0 

89.3 

2.3 

1028.7 

1053.3 

91.4 

1 . 9 

419 . 4 

442 . 3 

500.0 

17 . 2 

155697 

250000 

451.4 

(14.4) 

71536.1 

93922.6 

443 . 1 

(24) 

43727.4 

53650.5 

446.5 

1.3 

22772.4 

26333 .2 

459.1 

(  .3) 

8588.0 

11057.5 

DIFF 

(0-T) 


(44.3) 


(7.3) 


(3.1) 


(1.7) 


(1.1) 


(535.6) 


(14.9) 


3.5 


(2.3) 


(32.7) 


98 


(89.7) 


(24.6) 


(22.9) 


(94303) 


(9923.1) 
(3560 .8) 
(2469.5) 


Table  5.1:  Comparison  of  observed  to  theoretical  means  and 
variances . 
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difference  occurred  for  Expo (500)  and  the  simulated  variance 
was  60%  lower  than  the  theoretical  variance.  Both  differences 
are  mainly  due  to  the  fact  that  the  mean  times  between 
failures  are  quite  large  and  therefore  longer  simulation  runs 
will  be  necessary  before  we  could  expect  the  simulated  means 
and  variances  to  converge  to  the  theoretical  values. 

The  second  analysis  undertaken  was  to  determine  whether 
the  number  of  failures  cycled  over  the  years  of  the  engine 
life  cycle.  The  SIMAN  program  was  run  multiple  times  to 
simulate  the  potential  number  of  times  engines  failed  over 
each  year.  Appendix  G  provides  a  graphical  illustration  of 
the  output  from  this  process.  Table  5.2  adds  some  additional 
information  about  the  results. 

As  is  graphically  depicted  in  Appendix  G,  the  "curve"  of 
failures  has  the  greatest  "cyclic"  nature  when  the  exponential 
distribution  is  assumed.  When  the  Weibull  is  used,  even 
assuming  the  same  characteristic  life  parameter,  and  the  value 
of  beta  increases,  the  number  of  simulated  failures  in  each 
year  tends  to  smooth  out . 

This  can  best  be  interpreted  by  saying  that  for 
components  that  have  been  determined  to  be  "wearout" 
candidates,  and  after  a  Weibull  analysis  has  been  conducted  to 
determine  the  appropriate  parameters  for  the  Weibull 
distribution,  the  higher  the  beta  value,  the  more  likely  it  is 
that  there  will  be  a  more  constant  number  of  failures  per 
year.  In  the  case  of  an  assumption  of  a  high  characteristic 
life,  500  hours,  the  smoothing  effect  is  less  pronounced,  but 
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DISTRIBUTION 

AVG  NBR  OF 

COMPONENT 

FAILURES /YR 

- - - 

MAX  NUMBER  OF 

FAILURES  IN 

ANY  YEAR 

MIN  NUMBER  OF 

FAILURES  IN 

ANY  YEAR 

EXPO (20) 

120.7 

144 

102 

WEIB(20, 1.5) 

133.1 

144 

124 

WEIB (20,2) 

135.3 

144 

128 

WEIB  (20, 3) 

133.9 

139 

126 

WEIB (20,5) 

129 . 8 

133 

124 

EXPO  (50) 

47.2 

62 

35 

WEIB(50,1.5) 

52.4 

59 

48 

WEIB  (50, 2) 

53.4 

58 

50 

WEIB (50, 3) 

53 . 1 

59 

50 

WEIB (50, 5) 

51.7 

54 

48 

EXPO (100) 

22.0 

27 

11 

WEIBdOO,  1.5) 

25 . 1 

30 

20 

WEIB(100,2) 

25 . 8 

29 

22 

WEIBdOO,  3) 

25 . 1 

29 

22 

WEIBdOO,  5) 

25.3 

27 

21 

EXPO (500) 

5.3 

10 

0 

WEIB(500,1.5) 

5.2 

8 

2 

WEIB(500,2) 

5.3 

8 

3 

WEIB(500, 3) 

5.1 

7 

1 

WEIB(500,5) 

4.9 

8 

0 

Table  5.2:  STOimiiary  of  simulated  engine  failures. 
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is  still  present.  For  components  with  a  low  MTBF  (a  high 
failure  rate)  the  presence  of  the  smoothing  trend  is  quite 
pronounced . 

From  the  data  generated  by  the  SIMAN  simulation  program, 
ttie  author  concludes  that  the  Weibull  distribution  does  not 
show  an  increase  in  the  "cyclic"  nature  of  the  numbers  of 
failure  data,  but,  in  fact  shows  a  decrease  in  that 
phenomenon.  The  Weibull  distribution  for  beta  values  greater 
than  one  appears  to  have  a  "smoothing"  effect  on  the  number  of 
failures. 

The  third  part  of  the  analysis  was  to  determine  the 
nature  of  the  distribution  of  the  annual  number  of  failures 
when  the  time  between  failures  is  Weibull  distributed. 
Appendix  H  presents  histograms  which  can  be  used  to  identify 
some  characteristic  distributions.  The  data  for  these  plots 
was  developed  using  the  simulation  program,  and  they  represent 
the  number  of  failures  occurring  each  year  for  100  engine 
years .  The  means  and  variances  for  each  set  of  plots  are 
summarized  in  Table  5.3.  This  table  shows  the  mean  number  of 
failures  per  year  for  the  Weibull  decreases  initially  as  the 
beta  value  increases .  However  it  does  begin  to  increase  again 
once  the  beta  value  exceeds  approximately  2.2.  This 
phenomenon  occurs  for  the  same  reason  that  the  MTBF  decreased 
until  beta  reached  approximately  2.2  and  then  it  increased  as 
shown  in  Table  4.3.  The  variance  of  the  number  of  failures, 
however,  is  a  strongly  decreasing  function  of  beta. 

This  table  also  shows  that  the  means  and  variances  for 
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DISTRIBUTION 

MEAN  NBR  OF 

FAILURES /YEAR 

VARIANCE  OF 

FAILURES /YEAR 

EXPO (20) 

12 . 07 

10.22 

WEIB(20, 1.5) 

13.31 

5 . 02 

WEIB(20, 2) 

13.53 

2.87 

WEIB(20, 3) 

13.39 

1.83 

WEIB (20, 5) 

12 . 98 

.80 

EXPO (50) 

4.72 

4 . 14 

WEIB (50, 1.5) 

5.24 

2.34 

WEIB (50, 2) 

5.34 

-  1.09 

WEIB(50, 3) 

5.32 

.94 

WEIB (50, 5) 

5.17 

.46 

EXPO(IOO) 

2.2 

1.95 

WEIB(100,1.5) 

2.51 

1.32 

WEIB  (100, 2) 

2 . 58 

.83 

WEIB (100, 3) 

2.51 

.47 

WEIB (100, 5) 

2.53 

.29 

EXPO  (500) 

.53 

.49 

WEIB(500, 1.5) 

.52 

.35 

WEIB (500, 2) 

.53 

.31 

WEIB (500, 3) 

.51 

.25 

WEIB (500, 5) 

.49 

.27 

Table  5.3:  Mean  number  of  failures  per  year  and  variance 
of  nximber  of  failures  per  year. 
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the  exponentially  distributed  failures  are  close  in  value  but 
not  equal.  They  are  theoretically  equal  since  the  number  of 
failures  per  year  is  Poisson.  The  slight  differences  here  are 
consequences  of  a  limited  number  of  simulation  runs. 

On  the  first  page  in  Appendix  H,  for  those  graphs  with  an 
eta  value  of  20,  the  exponential  distribution  represents  an 
approximate  Poisson  distribution.  In  the  case  of  each  of  the 
"Weibull  plots"  they  do  not  look  like  a  Poisson  distribution 
or  even  a  discrete  form  of  a  Normal  distribution  when  beta  is 
less  than  two.  For  beta  values  greater  than  two  a  Normal 
distribution  using  a  continuity  correction  might  be  a 
reasonable  approximation. 

To  better  discern  the  nature  of  the  transition  from  the 
Poisson  for  an  eta  of  20  and  the  approximate  Normal  for  beta 
greater  than  two,  two  additional  plots  have  been  included  on 
the  second  page  of  Appendix  H.  These  plots  have  beta  values 
of  1.25  and  1.75,  and  an  eta  value  of  20.  The  mean  number  of 
failures  per  year  for  the  Weibull  (20 , 1 . 25)  is  12.87  with  a 
variance  of  7.12.  For  the  Weibull (20, 1 . 75)  the  mean  number  of 
failures  is  13.48  with  a  variance  of  5.24. 

The  histograms  for  eta  of  50,  100,  and  500  are  presented 
on  the  remainder  of  the  pages  in  Appendix  H  and  show  an 
increasingly  better  Poisson  distribution  for  the  exponential 
failure  times  than  was  the  case  for  the  eta  of  20.  The 
transition  to  the  Weibull  and  the  influence  of  beta  is  also 
increasingly  "smoother". 

The  fourth  aspect  of  this  analysis  is  to  determine 
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wh6th6ir  ths  Wsibull  Tna.y  hsvs  sny  impsct  on  financis.!  docision 
making.  That  will  be  addressed  in  the  next  section. 

C.  FINANCIAL  SPREADSHEET  ANALYSIS 

One  of  the  principal  goals  of  the  CIP  program  is  to 
maintain  an  engine  design  which  allows  the  maximum  aircraft 
availability  at  the  lowest  cost  to  the  government  [Ref.  10]  . 
Thus  the  question  remains,  "does  the  introduction  of  the 
Weibull  distribution  into  the  CEAMOD  program  more  accurately 
the  reflect  the  costs  and/or  potential  savings?"  By 
introducing  the  expected  failure  rates  per  year  into  the 
financial  spreadsheet  developed  by  the  author,  that  question 
can  be  answered. 

In  order  to  most  closely  reflect  the  General  Electric 
example,  shown  in  Appendix  C,  the  author  has  chosen  to  use  the 
exponential  MTBF  of  50  as  the  base,  and  compare  the  other 
three  values  (20,  100,  and  500)  to  that  base.  The  only 
exception  is  the  assumption  that  in  making  a  comparison  the 
designers  would  desire  to  move  from  a  lower  MTBF  to  a  higher 
MTBF.  Thus,  in  the  comparison  of  20  hours,  that  will  be  the 
"current"  and  50  will  be  the  "proposed".  In  the  other  two 
comparisons  50  is  the  "current"  and  either  100  or  500  is  the 
"proposed" . 

The  author  has  also  chosen  to  first  compare  similar 
distributions.  Thus  the  exponentials  will  be  compared  and 
the  Weibulls  with  like  beta  values  will  only  be  compared. 
After  that  a  comparison  will  be  made  between  the  exponential 
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and  Weibull  for  unlike  parameters.  In  that  comparison 
exponential  will  be  chosen  as  the  "current"  distribution  and 
one  of  several  Weibull  distributions  will  be  chosen  for  the 
"proposed" . 

The  highlighted  rows  in  Table  5 . 4  indicate  the  highest 
projected  saving  for  each  set  of  distribution  comparisons. 
Example  spreadsheets  from  which  this  data  was  extracted  are 
provided  in  Appendix  I.  The  savings  for  a  Weibull 
distribution  follow  a  similar  pattern  to  the  number  of 
failures  per  year,  increasing  as  beta  increases  to  two,  then 
slowly  decreasing  until  they  reach  the  level  achieved  by  the 
exponential . 

Based  upon  the  projected  costs /savings  table,  it  is  clear 
that  given  all  assuitptions  previously  discussed,  there  can  be 
a  significant  savings  differential  between  the  projections 
based  on  both  the  exponential  distribution  and  the  Weibull. 
However,  the  savings  between  Weibulls  is  larger  because  the 
difference  between  the  expected  number  of  failures  per  year 
for  the  current  configuration  and  the  proposed  configuration 
is  larger  with  the  Weibull  assxomption.  Thus,  one  can  expect 
a  greater  decrease  in  the  number  of  failures  per  year  with  a 
Weibull  assuit^tion  for  both  the  current  and  proposed 
configurations.  This  greater  decrease  translates  into  an 
increase  in  potential  savings.  For  example,  the  difference 
between  the  expected  annual  number  of  failures  for  Expo (20) 
and  Expo  (50)  is  7.35  (12.07  for  Expo  (20)  minus  4.72  for 
Expo(50)).  For  the  comparison  between  Weib(20,1.5)  and 
Weib(50,1.5)  the  difference  is  8.07  (13.31  for  Weib(20,1.5) 
minus  5.24  for  Weib ( 50 , 1 . 5 ) ) . 
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CURRENT 

PROPOSED 

PROJECTED 

FIRST  YEAR  OF 

DISTRIBUTION 

DISTRIBUTION 

SAVINGS 

POS  CASH  FLOW 

EXPO  (20) 

EXPO  (50) 

1 

WEIB  (20,1.5) 

WEIB  (50,1.5) 

$854,092 

1 

WEIB  (20,2) 


WEIB  (20,3) 


WEIB  (20,5) 


EXPO  (50) 


WEIB  (50,1.5) 


WEIB (50, 2) 


V«EIB(50,3) 


WEIB (50, 5) 


EXPO (50) 


WEIB (50, 1.5) 


WEIB(50,2) 


V1EIB(50,3) 


WEIB (50, 5) 


WEIB  (50,2) 


WEIB  (50,3) 


WEIB  (50,5) 


EXPO  (100) 


WEIB (100, 1.5) 


WEIB (100, 2) 


WEIB(100,3) 


WEIB (100, 5) 


EXPO (500) 


WEIB (500, 1.5) 


WEIB(500,2) 


WEIB(500,3) 


WEIB (500, 5) 


$867,797 


$854,092 


$824,397 


$220,221 


$244,205 


$247,632 


$253,342 


$233,926 


$410,953 


$471,485 


$481,764 


$481,764 


$466,917 


Table  5.4:  Projected  savings  from  changes  in  distribution 
parameters . 

Table  5.5  provides  a  comparison  of  the  projected 
costs/savings  when  changing  from  the  exponential  distribution 
to  the  Weibull.  This  might  correspond  to  a  situation  where  the 
old  component  had  a  constant  failure  rate  and  the  new  one  did 
not.  The  values  for  the  first  two  lines  are  negative  because 
the  comparison  is  between  distributions  with  the  same  eta 
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CURRENT 

DISTRIBUTION 

PROPOSED 

DISTRIBUTION 

PROJECTED 

(COSTS) /SAVINGS 

EXPO  (20) 

WEIB(20, 1. 5) 

(209,211) 

EXPO  (20) 

WEIB(20, 5) 

(171,522) 

EXPO  (20) 

WEIB(50, 1.5) 

712,470 

EXPO  (20) 

WEIB(50,2) 

701,049 

EXPO  (50) 

WEIBdOO,  1.5) 

184,816 

Table  5.5;  Expected  cost /savings  conparison  from  changing 
from  an  exponential  to  a  Weibull. 

value.  In  these  cases  the  average  number  of  failures  per  year 
increases,  thus  there  is  an  additional  cost  for  this  change. 
However,  the  cost  is  less  for  a  WEIB(20,5)  than  for  a 
WEIB(20,1.5)  because  the  niomber  of  expected  failures  per  year 
is  less  for  the  WEIB(20,5).  For  the  next  two  lines,  the 
comparison  between  an  EXPO (20)  and  a  WEIB(50,1.5)  and 
WEIB(50,2)  shows  a  savings,  but  not  as  great  a  savings  as  the 
change  shown  in  Table  5.4  when  the  change  was  made  with  like 
distributions.  This  is  due  to  the  increase  in  the  expected 
number  of  failures  per  year  between  distributions.  The 
difference  between  the  average  number  of  failures  for  a 
comparison  of  EXPO (20)  and  EXPO (50)  is  7.35.  However,  the 
difference  between  EXPO  (20)  and  Weibull ( 50 , 1 . 5 )  is  6.83  and 
for  Weibull (50, 2)  it  is  6.73.  This  decrease  in  the  difference 
between  the  current  and  proposed  distributions  accounts  for 
the  decrease  in  the  projected  savings.  The  savings  are 
primarily  achieved  from  the  shift  from  an  MTBF  of  20  hours  to 
one  approaching  50  hours  vice  a  shift  to  a  different 
distribution.  The  same  is  true  for  the  last  comparison  in 
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Table  5.5. 

Based  upon  this  analysis,  the  author  concludes  that  there 
is  a  significant  difference  in  costs  between  the  exponential 
and  Weibull  distributions.  In  addition,  when  applied  to 
Weibull  with  like  beta  values  but  differing  eta  values,  those 
savings  can  be  quite  large.  The  savings  for  a  Weibull 
distribution  follow  a  similar  pattern  to  the  number  of 
failures  per  year,  increasing  as  beta  increases  to  two,  then 
slowly  decreasing  until  they  reach  the  level  achieved  by  the 
exponential . 
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VI.  SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  SUMMARY 

The  key  questions  that  needed  to  be  answered  by  this 

thesis  research  effort  were: 

1.  Does  the  assumption  of  the  Weibull  distribution  mean 
that  there  will  be  "cyclic"  numbers  of  failures  over  the  years 
of  the  engine  life-cycle? 

2 .  What  does  the  probability  distribution  for  the  number 
of  failures  in  a  year  look  like  if  failure  times  are  Weibull? 

3 .  Is  there  any  significant  difference  in  the  net 
present  value  of  future  life  cycle  costs/savings  by  assuming 
the  Weibull  instead  of  the  exponential  distribution  for  the 
distribution  of  times  between  failures  for  an  aircraft  engine 
component? 

The  exponential  is  a  special  case  of  the  Weibull  and  it 
is  well  known  that  if  times  between  failures  are  exponentially 
distributed  then  the  number  of  failures  over  a  given  period  of 
time  are  Poisson  distributed.  However,  for  any  other  case  of 
the  Weibull  the  distribution  of  the  number  of  failures  over  a 
given  period  of  time  cannot  be  derived  analytically. 
Therefore,  the  use  of  simulation  is  needed  to  study  the  nature 
of  this  distribution.  As  a  consequence,  the  author  was  able 
to  generate  the  empirical  distributions  in  question  for  a 
range  of  Weibull  parameters.  And,  in  addition,  obtain  an 
estimate  of  the  mean  number  of  component  failures  per  year  for 
different  parameter  values  of  the  Weibull.  That  information 
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was  then  used  in  a  simple  version  of  the  CEA  model  to  compare 
the  results  derived  from  assuming  exponentially  distributed 
times  between  failures  with  the  results  from  assuming  Weibull 
distributed  times  between  failures. 


B.  CONCLUSIONS 

From  the  results  of  this  analysis  the  three  questions 
listed  above  can  be  answered.  The  answer  to  the  first 
question  regarding  "cyclic"  numbers  of  failures  is  that  the 
use  of  the  Weibull,  with  a  shape  factor  greater  than  one, 
leads  to  less  variation  in  the  annual  simulated  number  of 
failures  than  the  exponential.  As  the  beta  value  (shape 
factor)  increases  the  variance  in  the  total  number  of  failures 
is  reduced.  However,  the  Weibull  assumption  with  a  shape 
factor  greater  than  one  does  increase  the  expected  number  of 
failures  per  year  since  the  mean  time  between  failures  becomes 
smaller . 

With  the  data  analyzed  for  this  thesis,  the  distribution 
of  the  number  of  failures  per  year  for  the  Weibull 
distribution  appears  to  approximate  the  Normal  distribution 
with  beta  values  greater  than  two.  When  the  beta  value  is  one 
(the  exponential  distribution)  the  Poisson  distribution  is 
clearly  represented  and  can  be  theoretically  justified.  With 
beta  values  between  one  and  two,  the  distribution  appears  to 
be  transitioning  from  Poisson  to  a  Normal  approximation.  As 
the  beta  value  increases  the  skewness  of  the  distribution,  and 
the  variance  decrease  and  it  begins  to  approach  an  approximate 
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Normal  distribution  although  this  remains  to  be  statistically 
tested . 

The  question  regarding  differences  in  the  net  present 
value  of  the  life  cycle  costs  of  an  aircraft  engine  under  a 
proposed  component  change  also  produces  a  positive  result. 
If,  in  fact,  the  component  is  failing  in  a  manner  that 
reflects  a  Weibull  distribution  and,  there  is  a  proposal  to 
change  a  component  to  one  that  has  a  higher  MTBF,  the  use  of 
the  Weibull  may  show  a  greater  net  present  value  of  the 
savings  than  if  the  change  was  from  assuming  the  component 
failed  according  to  the  exponential  distribution.  Based  upon 
the  analysis  in  this  thesis  those  savings  may  be  significant. 


C .  RECOMMENDATIONS 

The  author  recommends  that  the  use  of  the  Weibull  in  the 
CEAMOD  program  continue  to  be  evaluated.  The  first  step  in 
this  process  is  to  determine  if  a  theoretical  distribution  can 
describe  the  number  of  failures  per  year.  A  Chi-square 
goodness  of  fit  test  of  potential  distributions  should  be 
made.  Next,  it  is  appropriate  to  determine  what  kinds  of 
components  actually  do  display  the  characteristics  of  the 
Weibull.  One  way  is  through  Weibull  analysis  as  shown  in 
Appendix  A.  If  a  component  has  its  time  between  failures 
described  by  a  Weibull  then  incorporation  into  the  model  is 
justified,  even  if  the  number  of  failures  distribution  cannot 
be  determined  precisely. 

In  order  to  incorporate  the  Weibull  into  the  CEAMOD,  the 
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progratu  designers  will  need  to  rewrite  portions  of  the  input 
and  calculation  sections  of  the  model.  A  failure  rate  can  not 
be  used  as  an  input  value  as  is  currently  done  for  the 
exponential.  The  expected  number  of  failures  in  a  given  time 
period  will  need  to  be  the  input  value  if  the  failures  are 
Weibull.  This  can  be  estimated  either  by  additional 
simulation  runs,  or  by  using  a  spreadsheet  generated  Weibull 
distribution  with  known  beta  and  eta  values  (determined 
through  Weibull  analysis)  as  described  in  Chapter  IV.  Once 
this  new  input  value  has  been  determined,  the  calculations  in 
the  financial  portion  of  the  model  can  be  modified  to 
accommodate  the  new  input  data. 

Finally,  this  thesis  has  also  considered  the  variance  of 
the  exponential  and  Weibull  distributions .  It  has  shown  that 
increased  beta  values  (greater  than  one)  result  in  reduced 
variance  of  the  number  of  failures  over  a  given  period  of 
time.  This  suggests  that  an  analysis  of  the  distribution  of 
life  cycle  costs  may  be  the  next  step  in  the  evolution  of  the 
CEA  model . 
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appendix  A;  WEIBULL  PAPER  AND  GRAPH 


This  first  page  of  this  appendix  contains  an  example  of 
blank  Weibull  graph  paper  taken  from  Reference  5.  This  paper  is 
used  to  plot  actual  failure  data  to  determine  whether  the  data 
characteristics  fit  the  Weibull  distribution. 

The  second  page  contains  a  sample  plot  of  three  points, 
which  is  full  explained  in  Chapter  III. 
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SAMPLE  WEIBULL  ANALYSIS 
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The  median  rank  table  assigns  the  median  value  to  each  item 
in  a  sample  size.  The  median  rank  is  used  for  plotting  against 
th  X  axis  of  Weibull  graph  paper.  The  highlighted  value  3,  is  used 
as  an  example  in  Chapter  111. 
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6.6 

16.2 

25.8 
35.5 

45.1 

54.8 
64  4 

74.1 
83.7 
93.3 


20 

3.4 
8.2 

13.1 
18.0 
22.9 

27.8 

32.7 

37.7 

42.6 
47.5 

52.4 
57.3 

62.2 
67.2 
72.1 
77.0 

81.9 

86.8 

91.7 

96.5 
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APPENDIX  C:  SAMPLE  CEAMOD  OUTPUT  REPORT 


This  appendix  contains  the  standard  three  page  output  report 
from  the  CEAMOD  program.  The  input  data  is  taken  from  Reference 
4,  and  is  the  example  data  provided  by  General  Electric. 
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TITLE 

ENGINE  MODEL  Fwwv-xx-yyy 
TASK/ECP 


F-zz 


CEA  VERSION  2  1  12/9/94 

Pg  1 


Task  Incorporation  Input 

1  0  Incorporation  Style  (1.2  or  3) 

1  =  Attrition 

2  =  Retrofit  at  1st  Opportunity 

3  =  Forced  Retrofit  Kits  /  Month 


2  0  Does  Kit  Cost  Replace  Normal  Maint  Material  Cost?  1=Yes  0=No 


3  0  Delta  Production  Cost 

4  0  Kit  Hardware  Cost  -  $  /  Engine 

5  0  Kit  Labor  Manhours  at  O&l 

6  0  Kit  Labor  Manhours  at  Depot 

7  0  Technical  Pubs  Cost  -  Total  $ 

8  0  TCTO  Cost- Total  $ 

9  0  Tooling/Support  Equipment  Cost-Total  $ 

10  0  Spare  Parts  Factor 

1 1  0  Scheduled  %  Events  being  Modified 

12  0  Unscheduled  %  Events  being  Modified 

1 3  0  Unscheduled  Event  Rate  allowing  Modification 

14  0  Production  Incorporation  Date 

15  0  Field  Incorporation  Date 


$10,000 

$15,000 

2 

20 

$500 

$1,500 

$500 

0% 

100% 

100% 

0020 

Year  - ^  1993 

Year  - ^  1993 


Scheduled  Input 

16  0  Scheduled  Maintenance  Interval  (TAC's) 

17  0  Calculated  Scheduled  Maintenance  Interval  Rate/1000  EFH 

18  0  Scheduled  Manhours  at  O  level 

19  0  Scheduled  %  Removed  at  O&l  level 

20  0  Scheduled  Manhours  to  Remove/Repiace  at  O  level 

21  0  Scheduled  Manhours  at  I  level 

22  0  Scheduled  %  at  O&l  requiring  Repair 

23  0  Scheduled  Repair  Cost  at  O&l  level 

24  0  Scheduled  %  Returned  to  Depot 

25  0  Scheduled  Manhours  at  Depot 

26  0  Scheduled  %  at  Depot  requiring  Repair 

27  0  Scheduled  Repair  Cost  at  Depot 

28  0  Scheduled  %  Scrapped 

29  0  Hardware  Cost  to  Scrap 

30  0  Scheduled  Engine  Test  Time 
Unscheduled  Input 

310  Unscheduled  Event  Rate/1000  EFH 

32  0  Unscheduled  Manhours  at  O  level 

33  0  Unscheduled  %  Removed  at  O&l  level 

34  0  Unscheduled  Manhours  to  Remove/Replace  at  O  level 

35  0  Unscheduled  Manhours  at  I  level 

36  0  Unscheduled  %  at  O&l  requiring  Repair 

37  0  Unscheduled  Repair  cost  at  O&l  level 

38  0  Unscheduled  %  Returned  to  Depot 

39  0  Unscheduled  Manhours  at  Depot 

40  0  Unscheduled  %  at  Depot  requiring  Repair 
410  Unscheduled  Repair  Cost  at  Depot 

42  0  Unscheduled  %  Scrapped 

43  0  Hardware  Cost  to  Scrap 

44  0  Unscheduled  Engine  Test  Time 

45  0  Unscheduled  Secondary  Damage  Costs 

46  0  Unscheduled  Incidental  Costs 

47  0  Number  of  P/N’s 
Optional  Input 

48  0  %  Improvement  in  Specific  Fuel  Consumption  from  Current  to  Proposed 

49  0  Aircraft  Loss  Rate  Improvement  /  1 ,000,000  EFH 


[  Standard  Inputs  ] 

Fiscal  Year  Dollars 

1991 

NPV  Rate 

10% 

Labor  Cost  /  Manhour  at  O&l 

$32  32 

Labor  Cost  /  Manhour  at  Depot 

$43.30 

Cost  to  introduce  new  P/N  -  $  /  PN 

$1,524 

Cost  to  Maintain  each  P/N  /  Year 

$250 

Fuel  Cost  /  Gallon 

$0  61 

Test  Fuel  -  Gallons  /  Hour 

150 

Flight  Fuel  -  Gallons  /  Hour 

150 

EFH  /  Year 

240 

TAC  /  EFH  Ratio 

300 

TOT  /  EFH  Ratio 

1  50 

Aircraft  Cost 

$0 

Month  _ ^  5 

Month  _ ^  8 

;  CURRENT  I 

PROPOSED 

]  3000  ; 

4000 

T 

I  1.000  I 

0  750 

1  00  I 

00 

]  100%  1 

100% 

I  10.0  ' 

100 

1  25,0  1 

25.0 

1  100%  1 

100% 

1  $500  1 

$500 

1  100%  1 

100% 

1  10,0  I 

100 

I  10%  1 

1% 

1  $25,000  1 

$20,000 

1  5%  I 

1% 

1  $62,500  I 

$50,000 

1  150  [ 

1.50 

I  0,020  [ 

0  002 

I  00  I 

00 

t  100%  J 

100% 

]  100  J 

100 

I  25.0  I 

25  0 

I  100%  j 

100% 

1  $500  1 

$500 

1  100%  I 

100% 

I  10,0  I 

100 

I  3%  I 

0% 

I  $1,250  1 

$1,000 

1  1%  1 

0% 

1  $62,500  1 

$5,000 

1  150  I 

1  50 

I  $100,000  1 

$100,000 

1  50  1 

$0 

1  ^  1 

4 

0% 

0  00 
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title 

engine  MODEL  Fwww-»<-yyy 
TASK/ECP 


STANDAFD  HtSTORY  FILE 


CEA  VERSION  2  1 


12/9/94 

Pg  2 


F-zz 


(N) 


Calendar 

Year 


(O)  (P) 

No  of  Available  Mod  Months 


(Q)  (R) 

Engine  Deitveries 


(S)  m 

Annual  Engine  Flight  Hours 

J  Average 
Fleet  1  per  Engine 


(U) 


Cumulative 

Engines 


(V) 

Attrition 


Cumulative 
Whole  Engine 


(W) 


Annual 

Whole  Engines 


title 

ENGINE  MODEL  Fwww-xx-yyy 
TASK/ECP 


F-zz 


CEA  VERSION  2  1 


1 2/9/94 
CEA  Guru 


ASSUMPTIONS 

a)  Incorporation  in  Production  engines  will  begin  in 

b)  Number  of  engines  produced  with  this  change  is 

c)  Number  of  spare  units  incorporating  this  change  is 

d)  Modification  of  operational  engines  can  begin  in 

e)  Incorporation  of  this  change  in  operational 
engines  will  be  accomplished  by  — > 

f)  Total  kits  installed  out  of  total 
engines  not  modified  in  production  is 

g)  Total  engines  lost  to  attrition  is 

h)  Total  engines  retired  unmodified  is 
1}  Estimated  yearly  flying  hours 


May 

1993 

550 

0 

Aug 

1993 

1st  Opportunity 

at 

Depot 

0 

of 

0 

50 

0 


240  EFH  /  Year 
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APPENDIX  D:  SIMAN  SIMULATION  PROGRAM 


The  first  page  of  this  appendix  contains  the  SIMAN  Model 
Frame  for  the  simulation  program  used  to  generate  the  Weibull 
data.  The  process  involved  in  the  program  is  full  explained  in 
Chapter  IV. 

The  second  page  contains  the  SIMAN  Experimental  Frame  used 
in  conjunction  with  the  Model  Frame  to  specify  the  specific 
parameters  needed  to  run  the  program.  The  specifics  of  the 
experimental  frame  are  fully  explained  in  Chapter  IV. 
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Begin, Yes , Yes ; 

/Model  for  the  Simulation  of  a  Weibull  Distribution  relative 
;to  the  Aircraft  Component  Improvement  Program 

;  This  is  the  Model  portion  of  the  SIMAN  program. 

;  Create  Entity  Module 

CREATE, 10, 0: , 1;  ' 

;  Creates  a  distribution  of  failures  of  engine  components 

;  ten  engines  at  time  zero,  max  of  one  batch 

ASSIGN:X (1) =X (1) +1; 

ASS IGN : ENGNUM=X ( 1 ) ; 

;  Assigns  a  sequential  number  1  to  10  to  each  engine 

Q1  ASSIGN:  Arrtime=TNOW; 

DELAY:WEIB(50,5) ; 

;  Delays  each  engine  for  the  operating  period,  end  of  delay 

;  is  the  failure. 

QUEUE, RepairQ, • 

SEIZE:Repairman,• 

DELAY :  0 ; 

RELEASE:  Repairman; 

;  These  steps  represent  the  repair  process,  delay  of  zero  to 

;  simulate  assumptions  of  CEAMOD  program 

TALLY :ENGNUM, INTERVAL (Ar r time) ; 

;  This  step  collects  the  data  on  the  time  of  failure 

DELAY: 0 :NEXT(Q1) ; 

END; 
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BEGIN; 

PROJECT,  Thesis,  Glenn  R.  Cook; 


ATTRIBUTES : 

RESOURCES : 
QUEUES : 
TALLIES : 


This  is  the  experimental  frame  of  the  SIMAN  program 
it  is  designed  to  allow  changes  in  the  data  for  different 
program  runs . 

Arrtime : 

ENGNUM; 

Repairman; 

RepairQ; 

1,  Time  for  Engine  1,  "Engl.dat": 

2,  Time  for  Engine  2,  "Eng2.dat": 

3,  Time  for  Engine  3,  "Eng3.dat": 

4,  Time  for  Engine  4,  "Eng4.dat": 

5,  Time  for  Engine  5,  "Eng5.dat": 

6,  Time  for  Engine  6,  "Eng6.dat": 

7,  Time  for  Engine  7,  "Eng7,dat": 

8,  Time  for  Engine  8,  "Eng8.dat": 

9,  Time  for  Engine  9,  "Eng9.dat": 

10,  Time  for  Engine  10,  "Engl0.dat"; 

REPLICATE,  10, 0, 240,N, Y, 0; 

This  indicates  that  the  program  will  run  once  for  2400 
units  of  time. 
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appendix  E:  SIMAN  OUTPUT  REPORTS 


This  appendix  contains  the  output  reports  generated  by  the 
SIMAN  program.  The  author  has  chosen  to  only  include  examples 
from  the  simulation  runs,  thus  this  appendix  only  has  the  reports 
from  the  runs  with  an  eta  value  of  20 .  The  reports  are  fully 
explained  in  Chapter  IV. 
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SIMAN  IV  -  License  #9050352 
Naval  Post-Graduate  School 


EXPO  (20) 


Project:  Thesis 

Analyst:  Glenn  R.Cook 

Run  execution  date  : 
Model  revision  date: 

11/27/1994 

11/27/1994 

Replication  ended  at  time 

:  2400.0 

TALLY  VARIABLES 


Identifier 

Average 

Variation 

Minimum 

Maximum 

Observations 

Time  for  Engine  1 

19.351 

.94969 

.20862 

83.182 

122 

Time  for  Engine  2 

17 . 520 

.  94349 

.35883 

73.424 

135 

Time  for  Engine  3 

19 . 741 

1.0464 

.12021 

103.66 

121 

Time  for  Engine  4 

21.066 

.95621 

.17511 

109.84 

113 

Time  for  Engine  5 

19 . 186 

.92286 

.09282 

96.620 

125 

Time  for  Engine  6 

19.296 

.83613 

.34302 

84.389 

124 

Time  for  Engine  7 

20.700 

.95878 

.27664 

116.35 

115 

Time  for  Engine  8 

20.483 

1 . 0106 

1 . 1141 

113 . 80 

116 

Time  for  Engine  9 

18 . 795 

.89282 

.10558 

90.237 

127 

Time  for  Engine  10 

21.799 

1.0095 

.07690 

108 . 14 

109 

Run  Time:  0  min(s)  2  sec(s) 
Simulation  run  complete. 


Mean  =  19.79  Hours 
Variance  =  355.7  Hours 
Fail/Eng/Year  =  12.07 
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SIMAN  IV  -  License  #9050352 
Naval  Post-Graduate  School 

Weibull (20,1.5) 

Project:  Thesis  Run  execution  date  :  11/27/1994 

Analyst!  Glenn  R. Cook  Model  revision  date:  11/27/1994 

Replication  ended  at  time  :  2400.0 

TALLY  VARIABLES 

Identifier  Average  Variation  Minimum  Maximum  Observations 


Time  for  Engine  1  17.710 
Time  for  Engine  2  17.271 
Time  for  Engine  3  16.077 
Time  for  Engine  4  19.058 
Time  for  Engine  5  16.727 
Time  for  Engine  6  19.110 
Time  for  Engine  7  19.280 
Time  for  Engine  8  17.493 
Time  for  Engine  9  19.460 
Time  for  Engine  10  17.619 

Run  Time:  0  min(s)  3  sec(s) 
Simulation  run  complete. 


Mean  =  17.98  Hours 
Varaince  =  143.0  Hours 
Fail/Eng/Year  =  13.31 


64641 

1.2111 

52.406 

135 

69992 

.55646 

62.322 

138 

65722 

.29150 

53.101 

149 

60095 

1.4608 

58.326 

125 

64935 

1.3782 

51.809 

143 

66742 

1.4340 

64.694 

125 

71869 

.  85501 

63.743 

121 

67730 

.59326 

55.313 

137 

64598 

1.1069 

59.900 

122 

69236 

.84961 

69.638 

136 
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SIMAN  IV  -  License  #9050352 
Naval  Post-Graduate  School 


Weibull  (20,2) 

Project:  Thesis  Run  execution  date  :  11/27/1994 

Analyst:  Glenn  R.Cook  Model  revision  date:  11/27/1994 

Replication  ended  at  time  :  2400.0 

TALLY  VARIABLES 

Identifier  Average  Variation  Minimum  Maximum  Observations 


Time  for  Engine  1  18.233  .52208  1.4531  50.979  131 

Time  for  Engine  2  17.019  .54696  1.4297  41.855  141 

Time  for  Engine  3  18.124  .48846  3.5033  39.402  132 

Time  for  Engine  4  17.884  .54376  1.2400  51.982  133 

Time  for  Engine  5  17.613  .52332  2.0428  46.506  135 

Time  for  Engine  6  16.724  .49560  2.3522  42.024  143 

Time  for  Engine  7  17.715  .48766  2.9432  45.533  135 

Time  for  Engine  8  18.128  .49125  2.7711  42.892  132 

Time  for  Engine  9  17.708  .54230  .83887  48.240  134 

Time  for  Engine  10  17.488  .50650  1.3625  45.068  137 

Run  Time:  0  min(s)  2  sec(s) 

Simulation  run  complete. 


Mean  =  17.66  Hours 
Variance  =  82.7  Hours 
Fail/Eng/Year  =  13.53 
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SIMAN  IV  -  License  #9050352 
Naval  Post-Graduate  School 


Weibull  (20,3) 

Run  execution  date  :  11/27/1994 

Model  revision  date:  11/27/1994 

Replication  ended  at  time  :  2400.0 


Project:  Thesis 

Analyst:  Glenn  R.Cook 


TALLY  VARIABLES 


Identifier 


Average  Variation  Minimum 


Maximum  Observations 


Time  for  Engine  1 
Time  for  Engine  2 
Time  for  Engine  3 
Time  for  Engine  4 
Time  for  Engine  5 
Time  for  Engine  6 
Time  for  Engine  7 
Time  for  Engine  8 
Time  for  Engine  9 
Time  for  Engine  10 


18 

.885 

.35075 

17 

.435 

.37529 

17 

.560 

.33070 

17  , 

.  641 

.38849 

17  . 

,803 

.33332 

17  , 

.500 

.37213 

18  , 

.  031 

.37175 

18  . 

094 

.36543 

17  . 

681 

.33839 

17  . 

921 

.33384 

2.4143  33.810 
3.3360  35.705 
5.2358  33.048 
4.1221  37.807 
6.5280  33.660 
3.1330  37.320 
3.4824  32.318 
5.6019  35.971 
4.8010  29.547 
4.1352  35.286 


127 

137 

136 

136 

134 
136 
133 

132 

135 

133 


Run  Time:  0  min(s)  3  sec(s) 
Simulation  run  complete. 


Mean  =  17.85  Hours 
Varaince  =  40.4  Hours 
Fail/Eng/Year  =  13.39 
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SIMAN  IV  -  License  #9050352 
Naval  Post-Graduate  School 


Weibull  (20 , 5 ) 

11/27/1994 
11/27/1994 

Replication  ended  at  time  :  2400.0 

TALLY  VARIABLES 


Project:  Thesis  Run  execution  date  : 

Analyst:  Glenn  R.Cook  Model  revision  date: 


Identifier 

Average 

Variation 

Minimum 

Maximum 

Observations 

Time  for  Engine  1 

18.216 

.21669 

9.0715 

27.425 

131 

Time  for  Engine  2 

18.339 

.26445 

6.9617 

29 . 079 

130 

Time  for  Engine  3 

18 .397 

.23844 

7.1913 

27 . 333 

130 

Time  for  Engine  4 

18 . 336 

.20049 

10.339 

26.597 

130 

Time  for  Engine  5 

18.592 

.20915 

6.5764 

27 . 137 

128 

Time  for  Engine  6 

18 . 693 

.21916 

8.6234 

28 . 126 

128 

Time  for  Engine  7 

18 . 137 

.22710 

6.8289 

27 . 794 

132 

Time  for  Engine  8 

18.587 

.21210 

9.3200 

27.573 

129 

Time  for  Engine  9 

18.436 

.20319 

7.7681 

27 . 680 

130 

Time  for  Engine  10 

18.332 

.22248 

7.7532 

28 . 317 

130 

Run  Time:  0  min{s)  3  sec(s) 
Simulation  run  complete. 


Mean  =  18.41  Hours 
Variance  =  16.6  Hours 
Fail/Eng/Year  =  12.98 
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APPENDIX  Ft  SIMAN  "DAT"  FILE  OUTPUT 


This  appendix  contains  the  output  of  a  SIMAN  "Dat"  file. 

The  first  coluran  is  the  cumulative  clock  time  from  the  first 
failure  to  the  total  run  of  2400  hours.  The  second  column 
represents  the  actual  times  to  failure  for  engine  1.  The  use  of 
this  data  is  fully  explained  in  Chapter  IV. 
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Time  for  Engine 
. 13679400E+03 
. 16538160E+03 
.30977810E+03 
.43940080E+03 
.55595120E+03 
. 62420570E+03 
.64869380E+03 
. 66302080E+03 
.69013670E+03 
.82652700E+03 
. 86262820E+03 
. 96482040E+03 
. 97945980E+03 
.10163370E+04 
. 11069960E+04 
. 12023170E+04 
. 15486620E+04 
. 16954530E+04 
. 17580750E+04 
. 19088880E+04 
.21682330E+04 
.22424920E+04 
.23284450E+04 
- . lOOOOOOOE+01 


1 

. 13679400E+03 
. 28587570E+02 
.  14439650E+03 
. 12962270E+03 
. 11655040E+03 
. 68254580E+02 
. 24488100E+02 
. 14326970E+02 
.27115910E+02 
. 13639030E+03 
.36101140E+02 
. 10219230E+03 
. 14639340E+02 
. 36877080E+02 
. 90659120E+02 
.95320680E+02 
.34634510E+03 
. 14679170E+03 
. 62621460E+02 
.  15081300E+03 
.25934500E+03 
.  74259520E+02 
. 85952880E+02 
. lOOOOOOOE+01 
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APPENDIX  G; 

GRAPHS  OF  ANNUAL  NUMBER  OF  FAILURES 


Year  Frequency 

1  4 

2  7 


EXPO  (500) 


I  FREQUENCY  FREQUENCY 


APPENDIX  H; 

FREQUENCY  OF  FAILURE  GRAPHS 


EXPO  (20),  FREQUENCY  OF  FAILURES 


0  2  4  6  8  10  12  1  4  16  18  20  22  I 

NUMBER  OF  FAILURES  PER  YEAR 


WEIBULL(20,1.5).  FREQUENCY  OF  FAILURES 


0  2  4  6  8  10  12  14  16  18  20  22 

NUMBER  OF  FAILURES  PER  YEAR 


WEIBULL  (20,2),  FREQUENCY  OF  FAILURES 


0  <N^<0®OWJ<0®S 

NUMBER  OF  FAILURES  PER  YEAR 


WEIBULL  (20,3),  FREQUENCY  OF  FAILURES 


WEIBULL  (20,5),  FREQUENCY  OF  FAILURES 
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number  of  failures  per  year 


number  of  failures  per  year 
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frequency 


EXPO  (50),  FREQUENCY  OF  FAILURES 


NUMBER  OF  FAILURES  PER  YEAR 


WEIBULL  (50,  J),  FREQUENCY  OF  FAILURES 


number  of  failures  per  year 


WEIBULL  (50.5),  FREQUENCY  OF  FAILURES 
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EXPO  (100),  FREQUENCY  OF  FAILURES 


2  3  4  5  6 

NUMBER  OF  FAILURES  PER  YEAR 


30 


2  3  4  5  6 

NUMBER  OF  FAILURES  PER  YEAR 


WEIBULL  (100,2).  FREQUENCY  OF  FAILURES 


NUMBER  OF  FAILURES  PER  YEAR 


WEIBULL  (100.3).  FREQUENCY  OF  FAILURES 


0  1  2  3  4  5  6  7 

NUMBER  OF  FAILURES  PER  YEAR 


WEIBULL  (100.5),  FREQUENCY  OF  FAILURES 


NUMBER  OF  FAILURES  PER  YEAR 
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EXPO  (500),  FREQUENCY  OF  FAILURES 


0  1  2  3  4  5 


NUMBER  OF  FAILURES  PER  YEAR 


0  1  2  3  4  5 

NUMBER  OF  FAILURES  PER  YEAR 


WEIBULL  (500,2),  FREQUENCY  OF  FAILURES 


0  1  2  3  4  5 

NUMBER  OF  FAILURES  PER  YEAR 


WEIBULL(500.3),  FREQUENCY  OF  FAILURES 


NUMBER  OF  FAILURES  PER  YEAR 


WEIBULL  (500.5),  FREQUENCY  OF  FAILURES 
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APPENDIX  I:  FINANCIAL  SPREADSHEET  OUTPUT 


This  appendix  contains  examples  of  the  output  from  the 
financial  spreadsheet  used  to  determine  potential  cost  savings. 
The  author  has  chosen  to  include  only  those  with  a  comparison  to 
the  mean  or  eta  value  of  20  hours.  The  output  from  these 
spreadsheets  is  fully  explained  in  Chapter  V. 
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